MÉMOIRE 
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Lu à l’Académie des Sciences, le 3 o avril 1817. 


L une des plus fécondes méthodes d’analyse que les géo¬ 
mètres aient imaginées, est sans contredit celle qui consiste 
à faire varier les constantes arbitraires, comprises dans les 
intégrales des équations différentielles, soit pour en déduire 
les solutions particulières, soit pour étendre ces intégrales 
à d’autres équations. Dans la théorie des perturbations des 
planètes , cette méthode se trouve indiquée par l’observation 
même, au moins en ce qui concerne les inégalités séculaires; 
car les constantes arbitraires sont alors les éléments du mou¬ 
vement elliptique, et d’après l’observation, ces éléments 
changent de siècle en siècle, excepté le grand axe qui reste 
invariable quand on fait abstraction des inégalités dépen¬ 
dantes de la configuration des planètes. Ces éléments sont au 
nombre de six : Euler a donné le premier des expressions 
différentielles de cinq d’entre eux, dans la pièce qui remporta 
le prix de notre ancienne Académie pour l’année 1766; mais 
ce n’est que vingt-cinq ans plus tard que la question a ete 
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traitée d’une manière directe et sous son véritable point de 
vue, par Lagrange qui a déduit les différentielles des six élé¬ 
ments elliptiques, du même principe auquel il avait ramené 
peu (Le temps auparavant les solutions particulières des équa¬ 
tions différentielles et l’intégration des équations linéaires 
qui contiennent un dernier terme indépendant de la variable 
principale. La matière semblait épuisée, lorsque cette impor¬ 
tante théorie reçut dans ces derniers temps une extension et 
un perfectionnement auxquels on était loin de s’attendre. La¬ 
grange etLaplace ont donné, en 1808, les différentielles des 
éléments elliptiques au moyen des différences partielles de 
la fonction perturbatrice prises par rapport aux éléments 
dont il s’agit et multipliées par des fonctions de ces éléments 
qui ne renferment pas le temps explicitement. Je n’ai pas 
besoin de rappeler l’avantage de cette heureuse transforma¬ 
tion, surtout pour le calcul des inégalités séculaires, des 
équations à longues périodes et de celles qui dépendent 
d’une cause spéciale, comme les inégalités du mouvement 
de la lune, dues à la forme de la terre. En suivant la direc¬ 
tion de son génie qui lui faisait saisir dans les résultats par¬ 
ticuliers , ce qu’il était utile de généraliser, Lagrange étendit 
ses recherches au mouvement d’un système quelconque de 
corps, sollicité par des forces que l’on 11’aurait pas considé¬ 
rées dans une première approximation ; et il parvint à des 
formules qui expriment les différences partielles d’une fonc¬ 
tion de ces forces, qu’on peut encore appeler la fonction 
perturbatrice, au moyen des différentielles des constantes 
primitives, multipliées par des fonctions de ces constantes. 
J’ai ensuite obtenu des formules inverses de celles-ci, qui 
donnent immédiatement les différentielles des constantes 
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arbitraires au moyen des différences partielles de la fonction 
perturbatrice, et dont les coefficients sont aussi des fonc¬ 
tions qui ne contiennent pas le temps d’une manière expli¬ 
cite. L’application que j’en ai faite successivement aux per¬ 
turbations du mouvement d’un point attiré vers un centre 
fixe suivant une loi quelconque, et à celles du mouvement 
de rotation d’un corps solide, m’a conduit à des formules 
identiques pour ces deux problèmes si différents l’un de 
l’autre ; de sorte que les constantes analogues dans les deux 
questions, ont la meme expression différentielle; résultat 
singulier que l’on peut considérer comme un théorème de 
mécanique. Ces deux questions comprennent toute l’astro¬ 
nomie : à la première se rapportent le mouvement des planè¬ 
tes , et celui de leurs satellites et des comètes ; à la seconde, le 
mouvement de la terre sur elle-même, troublé par l’action 
du soleil et de la lune, et la libration de notre satellite. Mais 
quoique ces deux problèmes de la translation et de la rota¬ 
tion des corps célestes, puissent ainsi dépendre d’équations 
semblables, les solutions qu’on en a données jusqu’à présent 
ne sont pas les mêmes; et c’est à faire disparaître cette diffé¬ 
rence que je me suis attaché dans ce nouveau Mémoire. Si 
l’on excepte une inégalité à longue période qui paraît affecter 
la longitude moyenne de la lune, mais dont l’existence 
n’est pas encore bien constatée, toutes les circonstances 
du mouvement des astres et de la terre que les observations 
ont fait connaître, les géomètres, et particulièrement l’auteur 
de la Mécanique céleste, en ont déterminé les iois et la 
cause d’après le principe de la gravitation universelle. Il ne 
reste guère maintenant qu’à simplifier les méthodes qu’ils 
ont employées; et c’est, en effet, les rendre plus simples et 
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les perfectionner que de les ramener autant qu’il est possible 
à l’uniformité. 

Dans le cas du mouvement des planètes autour du soleil, 
la petitesse des excentricités et des inclinaisons de leurs or¬ 
bites, permet de développer la fonction perturbatrice en 
une série de sinus des multiples de leurs moyens mouvements; 
or, on peut donner une forme semblable à cette fonction 
relative au mouvement de rotation de la terre, en observant 
que la terre tourne à très-peu près autour d’un de ses axes 
principaux, et considérant l’amplitude des oscillations des 
pôles de rotation a sa surface, comme une très-petite con¬ 
stante arbitraire dont, on aura à déterminer les variations 
dues aux forces perturbatrices. Cela étant, si l’on compare 
les six éléments arbitraires du mouvement de la terre autour 
de son centre de gravité, aux six éléments du mouvement 
elliptique, on aura d’une part cette amplitude et la longi¬ 
tude géographique de l’axe de rotation à une époque déter¬ 
minée, qui répondront à l’excentricité de l’orbite et à la 
longitude du périhélie; ensuite l’inclinaison de l’équateur 
et la longitude de son nœud sur l’écliptique, quantités ana¬ 
logues à l’inclinaison et à la longitude du nœud de 1 orbite; 
enfin, la vitesse angulaire de rotation et la longitude géo¬ 
graphique à l’origine du temps, d’une droite tracée dans le 
plan de l’équateur, qui remplaceront le moyen mouvement 
ou le grand axe dont il se déduit, et ce qu’on appelle, dans 
la théorie des planètes, la longitude moyenne de lepoque. 
C’est sous ce point de vue que je vais envisager la question 
qui fait l’objet de mon Mémoire. 
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' § 1 . 

Equations du mouvement de la terre autour de son centre 
dç gravité. 

(i) Le mouvement d’un corps solide autour d’un de ses 
points est complètement déterminé, lorsque l’on connaît 
à chaque instant la position de ses trois axes principaux pas- 
4 sant par ce point, celle de l’axe instantané de rotation par 
rapport à ces trois droites, et la vitesse de rotation du corps 
autour de cet axe. Nous distinguerons entre eux les axes 
principaux de la terre qui se coupent à son centre de gravité, 
par la grandeur des moments d’inertie qui s’y rapportent, et 
que nous représenterons par A,B,C, en supposant A le plus 
petit et C le plus grand. Cela étant, menons arbitrairement 
par le centre de gravité un plan fixe et une droite fixe, tracée 
dans ce plan ; au bout d’un temps t quelconque, soit 6 l’in¬ 
clinaison du plan des axes des deux moments A et B sur le 
plan fixe, <J> l’angle compris entre l’intersection de ces deux 
plans et la droite fixe, et <p l’angle que fait l’axe du plus 
petit moment A avec cette intersection ; ces deux derniers 
angles pouvant s’étendre depuis zéro jusqu’à quatre angles 
droits, et le premier seulement depuis zéro jusqu’à deux 
droits; et celui-ci étant toujours l’angle dièdre formé par les 
deux autres. Ces trois angles ô, et <p seront des fonctions 
de t; et quand leurs valeurs seront connues, la position des 
trois axes principaux, et, par conséquent, celle du sphéroïde, 
le seront aussi. 

Au bout du même temps t désignons par 8 la vitesse an- 
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gulaire de rotation, commune à tous les points de la terre, 
et par p, q, r, ses composantes relatives aux trois axes prin¬ 
cipaux ; de sorte qu’on ait 

8 2 =p' -hq 2 + r 2 , 

et que les rapports 

P L r 
8 ’ 8 » 8 

expriment les cosinus des angles que l’axe instantané de ro¬ 
tation fait avec les axes principaux du plus petit, du moyen 
et du plus grand moment d’inertie. Les trois quantités/?, q, r, 
seront aussi des fonctions de t; et l’on aura encre ces incon¬ 
nues et les précédentes, les équations 

rdt=dy — cos. Gcfy, 

q dt-=^s\w.§ cos. <%dty sin. <pc/0, (r) 

pdt= sin.Osin.' <p dty — cos. ©c/0, 

qui ne dépendent pas des forces appliquées au mobile, et 
résultent uniquement de ce que tous ses points forment un 
système invariable. Nous ne rappellerons pas l’analyse qui 
y conduit d’après cette considération ; nous ferons seulement 
remarquer que r étant la vitesse de rotation autour de l’axe 
du plus petit moment, rdt est l’angle décrit pendant l’instant 
dt par tous les points de la terre, parallèlement au plan de 
ses deux autres axes principaux; cet angle serait d^ si l’in¬ 
tersection de ce plan avec le plan fixe demeurait immobile; 
mais cette intersection décrivant un angle d ^ sur le plan 
fixe, dont la projection sur le plan mobile est cos. il 

en résulte, en ayant égard au sens de l’accroissement des 
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angles t|* et <p, que rdt est l’excès de cLy sur cos. 6^^; ce qui 
donne la première équation (i). 

Soient encore M, M', M", les sommes des moments des 
forces qui agissent sur tous les points de la terre, rapportés 
respectivement aux axes du plus grand, du moyen et du plus 
petit moment d’inertie; nous aurons ces trois autres équa¬ 
tions : 

Cdr+ (B— A)pqdt—Mdt, \ 

Bdq-h(A — C)prdt=M'dt , J (2) 

Adp -4- (C — B) qrdt — Wdt, ) 

qu’il faudra joindre aux précédentes. 

Le problème consistera donc à intégrer ces six équations 
différentielles du premier ordre, pour en déduire les valeurs 
des six inconnues p, q, r, ô, <|/ et 9 en fonctions de t et de 
pareil nombre de constantes arbitraires; ce qui sera possible 
par des approximations successives, fondées sur la petitesse 
des forces perturbatrices de la terre dans son mouvement 
de rotation. Nous n’irons pas au-delà, dans ce Mémoire, des 
quantités du second ordre par rapport à ces forces. En les 
négligeant tout-à-fait dans une première approximation, les 
équations (2) se réduiront à 

Cdr -+- (B— A)pqdt=o, 

BJ<7-t-(A— C)prdt—o , ( 3 ) 

A dp -h (C — B) <7r^* = o. 

(2) On sait intégrer les équations (1) et ( 3 ) sous forme 
finie, au moyen de deux intégrales définies; mais dans la 
question qui nous occupe, il vaudra mieux employer des 
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intégrales approchées de ces équations, qui nous permet¬ 
tront d’exprimer explicitement les inconnues en fonctions 
du temps, ce qui ne serait pas possible d’après leurs inté¬ 
grales exactes. Cette approximation sera fondée sur ce fait 
que les déplacements des pôles de rotation de la terre à sa 
surface sont très-petits, puisque les observations n’ont pu 
jusqu’à présent les rendre sensibles. Cela ne peut avoir lieu 
qu’autant que l’axe de rotation s’écarte très-peu de l’un des 
axes principaux, qui doit être celui du plus grand moment 
d’inertie, à cause que la terré est un sphéroïde aplati sur 
ses pôles. Il en résulte que les variables p et q sont très- 
petites par rapport à r ; en négligeant leur produit, la pre¬ 
mière équation ( 3 ) donne d’abord 

dr~o, r—n; 

n étant une constante arbitraire, qui représentera à- très-peu 
près la vitesse de rotation de la terre. Nous la supposerons 
positive, ce qui est permis, et revient à dire que l’angle 9 
est compté dans le sens de la rotation, l’inclinaison 6 étant 
moindre qu’un angle droit : l’angle ^ se compte en sens con¬ 
traire à partir de la ligne fixe qui en est l’origine. Si l’on 
substitue n à la place de r dans les deux autres équations (3), 

11 vient 

B dq — (C — A)npdt—o , 

A dp + (C — lï)nq dt=o. 

En intégrant celles-ci, faisant pour abréger, 



et désignant par e et c les deux constantes arbitraires, nous 
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aurons 

p=aes\n.ab(nt c), 

<? =—b ecos.ab(nt -b c). 

Les quantités a et b étant réelles par hypothèse, il suffira et 
il sera nécessaire que la constante e soit très-petite pour que 
chacune des variables/; et q le soitausri, comme nous lavons 
supposé. 

Je substitue ces valeurs de p et q dans la première équa¬ 
tion ( 3 ), d’où je déduis ensuite 
(B — A)e a 

r— n — cos. 2 ab(nt-bc), 

Cette seconde valeur de r étant substituée à son tour dans 
les deux dernières équations ( 3 ), on en déduira de nouvelles 
valeurs de p et q, plus approchées que les précédentes, et 
en continuant ainsi, il est évident qu’on obtiendra des va¬ 
leurs de p, q y r, ordonnées suivant les puissances de e. 
On voit aussi que les valeurs de p et q ns contiendront que 
des puissances impaires, et celle de r des puissances paires; 
si donc nous négligeons le cube et les puissances supérieures 
de e y il faudra s’arrêter aux valeurs de p, q y r, que nous 
venons d’écrire. 

( 3 ) Après les avoir substituées dans les équations (i), on 
en déduit 

/g_ e * 

cos. iab{nt -f- c)dt -f- cos.8 

dü = — [£sin. <pcos.ab(nt -h c) -b a cos. vs\n.ab(nt cj\edt y 
sin.Gfi?<j> = — [&cos.çcos.<z6(rc t-bc )—asin.çsin.a&(/2f + c)\edt. 

En négligeant la quantité e, et désignant par *, g et / trois 
constantes arbitraires, on aura d’abord 
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0=^, tyz=g> <p =z7lt-t~l. 

Si l’on tient compte de la première puissance de e, que l’on 
ait égard aux valeurs de a et &, et que l’on fasse, pour abréger, 

~ cos. (n t -h V) cos. a b (n t H- c)— ^ sin. (n t + /) sin. a b (n t + c)=k, 

H sin. (n t -+- 1 ) cos. a b {n t + c) H- ^ cos. (nt + l) sin. a b (ji t -h c) = k', 

on trouvera, dans cette seconde approximation, 


0 =>. 


k e 
n 


= fl t H- l 


*t = ë 

k‘ e cos. X 
«sin.^ 


k' e 

n sinA ’ 


Soient ensuite £0, et $9, les parties de 6, 4» et <p qui dé¬ 
pendent du carré de e; en faisant, encore pour abréger, 


g~-|[B(C—A) A(C— B)][i — cos.o.(nt + l)cos.2.ab (nt + cj\ 
— C(B — A) [cos. z{nt + l )—cos. 2.ab(nt -+-c)\ 

2 AB ab sin. 2 (nt q- 1 ) sin. 2 a b (nt -t- c) j = h, 


j[B(C — A) + A(G — B)] sin. 2 (nt H- l)cos.2ab(nt-+ c ) 
G (B — A) sin. 2 (n t -t- V) 

2 AB<3#cos.2(rc£ + /)sin. 2ab(nt + c)\ =h', 


on trouvera facilement 

. he'cos.'k 


/i’sin.X * 

_(aC—-A—B) e*t 

* 4^C 


_a h' e* cos. 1 

«’sin.’X 1 
h' g 3 (i +cos. 3 X) 
rc a sin. a >. 
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Il est évident quen continuant ces approximations suc¬ 
cessives, on formera des valeurs de 6, et 9 en séries ordon¬ 
nées suivant les puissances de e : nous nous arrêterons, comme 
plus haut, au carré inclusivement; mais afin que n soit la 
seule constante multipliée par t, nous comprendrons le pre¬ 
mier terme de &<p dans le terme ni de <p, ou, autrement dit, 

nous mettrons n — à la place de n; ce qui mo- 

difîera les valeurs de r et <p, et ne changera rien à celles de 
p t q, 0 et ^, puisqu’on néglige les puissances de c supérieures 
à la seconde. De cette manière, la solution complète des six 
équations (1) et ( 3 ) sera exprimée par ces formules : 


p — aesiu. a b (n t 4- c), 
q = — b e cos. a b{nt -h c), 


„_rc__ (2C-A-B) g 


1 —- 




4«c 
he^cos.y. 

+* J—: s- , 

rrsin.A 7 

k'e 2 A'e’cos. A 


(B—A)e a 7 . 

■ COS. + 


<p = nt-h l 


«sin.A « a sin. a A 7 

k'ec os. A . h' e a (i -h cos. a A) 


«sin.A 


sin. a X 


(4) 


qui comprennent les six constantes arbitraires n, e, c, L 
( 4 ) Maintenant pour avoir égard aux forces perturbatrices, 
et étendre cette solution aux équations (i) et (2) qui les con¬ 
tiennent, nous emploîrons la méthode connue de la varia¬ 
tion des constantes arbitraires. Les formules qui expriment 
leurs variations, renferment deux fonctions que nous repré¬ 
senterons par ü. et T, dont la première prise avec un signe 
contraire, sera l’intégrale de la somme des force perturba- 
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trices, multipliées chacune par l’élément de sa direction, et 
la seconde désignera la demi-somme des forces vives de tous 
les points du sphéroïde. On y regarde comme une fonc¬ 
tion donnée de t et des constantes arbitraires, et T comme 
une fonction aussi donnée des variables indépendantes qui 
fixent la position du mobile à chaque instant, et de leurs 
différentielles par rapport à t. Dans le cas actuel, ces varia¬ 
bles sont les trois angles 9 , <[> et 9 ; en faisant donc 


dt ~ 


. .dû do , 

9 ’ ^=+ ’ 


on supposera qu’on ait exprimé T sous la forme : 

T = fonct.(0, ^,9,0',<|/, 9'). 

Cela étant, faisons 


dT y dT _ d T_ 

dÿ dÿ — *' 

de sorte que £, 71 et w soient, ainsi que 0 et 9 des fonc¬ 
tions données de t et des six constantes n, c, e, \ , l. Dé¬ 

signons pardeux de ces constantes; et soit 

, f _ d^_d\ _ d^ d^ d§_ djn _ d§_ d’D 

V’f ] ~~dfdf df’ df^~ df df df df 

d<% d(à d<p d là 

+ Tjdf~drdf'’ 1 ' 

notation d’après laquelle on aura 

(//)-<>• 


Les formules d’où dépendent les différentielles des constantes 
arbitraires devenues variables, seront de la forme *. 
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-2j.dt={f,n)dn + {fe)de + (fc)dc + (f,\)d\ 

+ (fig)dg + (fil)dl. ( 6 ) 


On y mettra successivement n,e y c,\,g, l, à la place de f, 
ce qui donnera six équations linéaires, d’où l’on tirera les 
valeurs des six différentielles de ces constantes. Il existe 
d’autres formules inverses de celle-ci, qui donnent immé¬ 
diatement les valeurs de ces différentielles au moyen des 
différences partielles de Cl par rapport à ces mêmes con¬ 
stantes, mais pour lesquelles le calcul des coefficients serait 
moins facile que celui des quantités que nous désignons par 
C ff)- On sait, au reste, que ces quantités doivent être des 
constantes absolues ou des fonctions d’une ou plusieurs des 
six arbitraires, qui ne renferment pas le temps explicitement ; 
et c’est, en effet, ce qu’on va vérifier par le calcul de leurs 
valeurs. 

( 5 ) A un instant quelconque, soit r, la distance d’un point 
quelconque du sphéroïde à l’axe instantané de rotation , et 
z, ses trois coordonnées parallèles aux axes principaux 
du plus petit, du moyen et du plus grand moment d’inertie. 
A cause que les cosinus des angles que fait l’axe de rotation 


avec ces trois droites, sont respectivement g , |, g » on aura 


r I ' = x l 1 +J, 1 -h z* 



tZi 

« ' 


r*,y 

TJ • 


D’ailleurs en désignant par dm l’élément différentiel de 
la masse, et intégrant dans toute l’étendue du sphéroïde, 
nous aurons aussi à cause des axes principaux : 
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Jipc* +y, 2 )dm==C, J{z; + x I ')dm=B, JiyS-h z, 2 )dm= A; 

au moyen de quoi, nous conclurons de 1 équation prece¬ 
dente : 

/ 'Ïr^dm—Ap 2 -h B q 2 -t- C r 2 ; 

or, cette intégrale exprime évidemment la somme des forces 
vives de tous les points du sphéroïde; nous aurons donc 
pour la demi-somme 

T=±Ap 2 +' 7 Bq 2 + '-Cr 2 . 

D’après les équations (i), on a 

yp=z<j/sin.8sin.<p—G'cos.ç, 

2 = <|/sin.6cos.9 h- G'sin.9, 
r=( p' — <|<cos.&'; 

si donc on difïérentie T par rapport à <|/, 8 ' et 9 ', il en ré¬ 
sultera 

£ = A/>sin.Gsin.<p + B<7sin. 8cos. 9—C/’cos.6, 
7)=B^sin. 9 — A p cos. 9, 

«=C r ; 

et en substituant dans ces formules, à la place de p, q , r, 8, 
et 9, leurs valeurs données par les équations (4), on aura 
celles de » et w dont on a besoin, c’est-à-dire en fonctions 
de t, n , e, c> \,g et l. 

La substitution faite, et en négligeant toujours les puis¬ 
sances de e supérieures au carré, on trouve 

S = — Crccos.x + — B )l ‘- cos ^ 

2Ü ? 
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V) 


— C ek' H- 


2 CA'e 2 cos.>. 
/i sin.X 


^ (2C — A—B)e a (B — A)c 2 , . \ 

a = — v -— -*——cos. 2.ab(nt-\-c). 


On peut remarquer que la quantité £ est elle-même une con¬ 
stante arbitraire; ce qui tient à ce que \dt exprime la somme 
des projections sur le plan perpendiculaire à l’axe du plus 
grand moment d’inertie C, des aires décrites pendant l’instant 
d t par toutes les molécules de sphéroïde, multipliées par leurs 
masses respectives. 

(6) Il faut donc mettre ces valeurs de £, n et to, et celles 
de ^ , 0 et <p dans l’équation ( 5 ) pour en déduire les valeurs 
d ) relatives aux six constantes prises deux à deux. 
Mais quoique nous ayons tenu compte jusqu’ici des quantités 
du second ordre par rapport à e, les différences partielles 
de £, yi, etc. relatives à cette quantité, et par conséquent les 
valeurs de (n , e), (c, e), etc. qui en dépendent, ne seront 
exactes qu’aux quantités près de cet ordra exclusivement ; 
c’est pourquoi nous nous bornerons à ce degré d’approxima¬ 
tion pour toutes les valeurs de (/•,/'), c’est-à-dire, que nous 
rejeterons dans ces valeurs les termes dépendants du carré et 
des puissances supérieures dee; ce qui en rendra le calcul 
plus facile. 

On trouve alors neuf de ces quantités égales a zéro, savoir, 

("• é) = o , (re, c) = o, {n , \) = o, (e, X) =o, (c, \)= o,(c,g) = o, 
(c , l) = o, (>, /)= o, (g, l) = o ; 

et pour les autres, il vient 

{n,g) = Ccos.X, 

(/!,/)= C, 
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{e,c)-= — 


(C — A)(C—B)e 


/. -V_ (C-A)(C-B)«o S .X 

-^C”-’ 

(X,g‘) =— Cn sin.X. 

(7) Au moyen de ces résultats, 1 équation (6) donne : 


^f=Cc6s.X dg— Cdl, 

^dt= ^=^§^ (edl—edc—ecos.y.dg), 

i te £zi=!) erfe) 

de n G 


r dt= — C n sin. X dg, 


i n ^ j (C — A) (C—B)cos.X j ^ 7 

-dt= —Ccos .\dn -4---^-7;— L - ede -+- Crcsin.Wx, 

«Li 1 


-d tz=Gdn- 


(C-A)(C-B), 


d’où l’on tire réciproquement 

j / dci dn\ dt 
dn = {lTc + 7ï) c> 


eue— de (G— A)(C— B) ’ 
j dn Cndt daedt 

eac — ^ ( C __ A ) (C —B) Ta ~C 

Jmt dn dt dn cos Xdt 
^ dg C n sin. X dl Cn sin. X ’ 

1 __ dn dt 

g d\ C n sin.X ’ 

,1 rfnCrcsin.X dçidt 

dX cos Xdt dn C 
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D’après les valeurs de 0, et 9 que l’on devra mettre dans 
cette quantité sera une fonction de n y nt-\- c et nt+l; 

si donc on désigne pour un moment par la partie de 

la différence partielle, relative aux n non comprises dans 
nt + c ou nt-\- 1 , sa valeur entière sera 


do _/a?n\ 

du V dnj 


£n £a . 


et en vertu des équations précédentes, on aura 


dl 4 - tdn = - 


d 0 Qitidt 
~ ~de (C- A) (C — B) 
do cos.\dt / dci\dt 
~ dl C/isin.X \dn) G 


J J v do Cndt /'«UN edt 

e(dc + tdn) — — ~j- e ^ c _ ^ (C —B) “c" ’ 

J. ^ j do cos.\dt / dci\ dt 

dl -t- 1dn = — ^ CÏÏÏÎÏÏ3 — \25 ) C • 

Il en résulte que si l’on remplace c et respectivement par 
c—Jtdn et l — f*tdn, et si l’on fait attention qu’on a iden¬ 
tiquement 

nt — Jtdn—Jndt, 

la quantité a se changera en une fonction (kn,c+ Jndt 

/4- j'ndt , et les différentielles des nouvelles quantités c et l 

seront toujours données par la troisième et lp, sixième équa¬ 
tion (7), dans lesquelles les différences partielles relatives à 
n se rapporteront seulement à cette variable en dehors du 

signe J', en sorte que si nous faisons 

J'n dt — pi, 


elles seront prises en regardant la quantité comme con- 
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stante , ce qui en simplifiera l’expression. En même temps 
nous aurons 

dg ^ dg dg . 

de dl d\ l 7 

et la première équation (7) deviendra 

j dg dt zo\ 

( g ) 

On pourra aussi transformer plusieurs des quantités n , e, c, 
X,g-, l en d’autres, et exprimer les différentielles de celles-ci 
au moyen de différences partielles de Cl qui s’y rapportent. 
Soit, par exemple , 

esin .abc-=zf i ecos. abc'—f'\ 

on en conclura 


df=s\n.abcde -+- abecos.abcdc, 
df = eos .abcde — abesxn.ab c de > 

d g _ f dg f'dg 

de e df + e df' ’ 
d g » d g j d g 

dï= ah Sdf- ah fdf-’ 

et d’après la deuxième et la troisième équation (7) , on aura 
pour les remplacer : 

j,. dg C ndt dg abf'dt \ 

df'XKab ~df C 5 | 

1 w_ dg Çindt dg abfdt ( 9 ) 

~-~dfÂB^b + d7i~~(r~ ’ J 

en observant que 

(C—A) (C—B) = ABæ’ b\ 
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Ces différentes formules, jointes aux équations ( 4 ), ren¬ 
ferment la solution demandée des équations (i)et (2). Nous 
allons développer les conséquences qui en résultent, soit 
relativement aux pôles et à la vitesse de rotation de la terre, 
soit par rapport au mouvement de ses axes principaux autour 
de son centre de gravité. Ces formules conviendraient aussi 
au mouvement de rotation de la lune, troublé par l’action 
de la terre; mais on n’en déduirait par les mêmes consé¬ 
quences que pour le sphéroïde terrestre, parce que dans le 
cas de la lune, il existe un rapport commensurable, le rap¬ 
port d’égalité, entre le moyen mouvement de rotation et le 
moyen mouvement du centre de la force perturbatrice. 

§ II. 

Permanence de la vitesse et des pôles de rotation. 

(8) Les dimensions de la terre étant petites relativement aux 
distances des astres qui agissent sur ses différents points, on 
pourra toujours réduire H en une série convergente, ordonnée 
suivant les puissances négatives de ces distances. Les termes 
de cette série renfermeront des sinus ou cosinus des mul¬ 
tiples de ^ et de a bp, provenant des valeurs de 0, + et <p, et 
ils varieront en outre en raison du mouvement des astres 
auxquels on aura égard. Ils seront tous très-petits; et parmi 
les termes correspondants des valeurs de dn , df, etc., on 
devra rejeter ceux dont la période sera d’un jour ou d’un 
sous-multiple du jour, et ne conserver que les termes con¬ 
stants. s’il en existait, ou bien les termes à longues périodes 
qui acquerront par l’intégration de petits diviseurs, capa- 

4 
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blés de les rendre sensibles dans les valeurs de n,f 9 etc. Gela 
posé, occupons-nous d’abord de la formule (8). 

Si l’on ne tient compte que de la première puissance des 
forces perturbatrices, et si l’on regarde conséquemment dans 
cette formule, n,f, etc., comme des constantes, la différentia¬ 
tion relative à ^ n’y laissera subsister que des termes périodi¬ 
ques dépendants de nt, ou d’autres termes dépendants de 
abntet affectés du facteur a b, de sorte que les uns et les autres 
resteront insensibles après l’intégration. A la vérité, la valeur 
numérique de ab est telle que abnt ne surpasse le moyen 
mouvement du soleil que d’environ un septième; il en résulte 
que les termes qui auront pour argument abnt moins ce 
moyen mouvement seront sept fois plus grands dans l’inté¬ 
grale de la formule (8) que dans l’expression de ~ ; mais 

on verra bientôt quils sont si petits dans le devéloppement 
de H, que cette circonstance particulière ne peut les rendre 
aucunement appréciables dans la valeur de n, ni même dans 
celle de p qui s’en déduit par une seconde intégration. 

Nous pouvons donc, sans crainte d’erreur, considérer 
comme constante la valeur de n qui résulte de l’équation (8). 
S il ne s agissait que de cette quantité, il ne serait pas nécessaire 
de considérer les termes de iî, d’un ordre supérieur au pre¬ 
mier par rapport aux forces perturbatrices; mais il importe 
de faire voir que les termes du second ordre ne sauraient 
non plus devenir sensibles dans la valeur de quoiqu’ils 
y subissent deux intégrations successives. 

(9) Pour le prouver, désignons par 8 f, etc. les 

parties de j*, ».,/*, etc., qui sont du même ordre que ces 
forces, et par et ï'/t les parties de ^ et n qui.soat du 
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second ordre et qu’il s’agit de considérer. D’après l’équation 
(8), nous aurons 


C d.iïu _ 


d'il , 


dt 


d ^ 
d'n 


d'il 

d\j.dn 


% n 


d'n 

d\Ldf 

d 


V 


d^df d[r.d\^ d\Ldg 


ci * d" 1 n N 

-*e + dïdi* 


On regardera dans ces différences partielles secondes, les 
quantités n comme constantes, et l’on y fera 
pr —nt; on traitera de même les seconds membres des équa¬ 
tions (7), (8) et (9) pour en déduire par l’intégration, les 
valeurs de In, If, etc. ; de cette manière, on aura 


r d.8*n i d'n f('dn . 1 'd'n fdn j A 

0 dt —c dpj) dp ar ^cd^dTj 7^ at 
. c n ( d'n fdcij. d'n Cdn 
-*-ABab\d l Ld/'J df at d^dfJTf at ) 

ab_ r r d'n _ w d'n _1 d'n \ rdn , 

C \y d\j. df r ^ d\^df a b d\j.dï) J dn * 

C n sm.X \dÿ.dkj dg d\udgj d\ J 




C/zsin.A \d\j.d\J dl 


D’ailleurs, on a identiquement 

~ t. f ri dci fdn \ dn dn dn 
ab V If—fdJ-'J^-Tc' Tc+-d7 = 

le coefficient de f~dt se réduira donc à 
J dn 


dn 

dp’’ 


d'n 


C d [l 7 

et la formule précédente pourra s’écrire ainsi : 


4- 
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J /• /*Ja 

(io) 


d.Ù'n 

dt 


=<?$//!? **+ R ; 


oii i’on représente par R une quantité composée de parties 
qui seront toutes de la forme : 

.(^Lf\ p dt—lJLp P ’dt) ; 


Y W| d v-'_ 


Y étant un coefficient constant, P et P' désignant deux fonc¬ 
tions de t, pour lesquelles on prendra successivement 


d n dtt dfi dit 
v et dn » ~df et dj~' ’ 


d J/. 


dçi dct_ d Cl dçi 

~d~g et d\ ’ HT et ~d\' 


Cette formula ( io) se compose ainsi de deux parties dis¬ 
tinctes que nous allons examiner l’une après l’autre. 

(io) Si l’on fait abstraction des termes de la première partie 
dont les périodes seront d’un jour ou d’un sous-multiple 
du jour, pour ne conserver que ceux dans lesquels les mul¬ 
tiples de nt auront disparu, il est évident qu’il faudra em¬ 
ployer en dehors et sous le signe JJ, des termes de dé¬ 
pendants d’un même multiple de nt. Soit donc 


Q. = YlcQS.(int-bj abnt -—a t —ê) 

-t H'cos f abnt — a 't -—£'); 

i,j,f étant des nombres entiers ou zéro, H,H', a, a',ê, 
des constantes données, et a£, a* t, provenant des mouve¬ 
ments des astres qui agissent sur la terre. En négligeant, 
comme nous le disons, les termes à courte période, nous 
aurons 
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■HH'(i+/ai)(i+/'ai)r ffî 


i-\rjab 


(int+jabnt—ct)' 


(i.n -f -j 


(a- 


fq fil ;^‘>'iît>fi soi 


et d’après la première partie de la formule (io), on aura en 
même temps 


= _ H JL r 

4 * (i hj—J') ab +l (« — *') 

— T j cos. [(j‘~j)abnt -4-'(à — a") ï-h '6'—g]; 


ou l’on a fait pour abréger 

f Vi — - • [ ï ' •' . = i| - C* 

(i -h j a b) (i+j’ab) [ r—] --7 Vv1 = ï ^ 


( i -t -j a £)(*+/' a £) [i i + ( / ■*-/) a b] [a i -t- (j -t -/) ab—*(<x-\- a'JJ 


I 


-V?- 

Q+jab — ^ • ; TTti'î >DG-VI 

' fi qj'rn^irr: l • i • > 

Or, si i n’est pas nul, on aura à peu près I —a et 1=^4 


à cause que ab, ^ et ^ sont de petites fractions; si l’on a 

i= o, ces quantités I et I f , aussi bien que la fraction ijui 

multiplie I, ne pourront devenir très-grandes, a moins que 

>r 'duffil U . ‘jftori 

parmi les valeurs numériques de jab — -, il y en ait une ou 

plusieurs qui soient de très-petites fractions dey ab; ce qui n’a 
pas lieu , puisque le cas le plus fayprable serait celui où l’on 

aurait /= f et ^ égal au moyen mouvement du soleil, cas 
èyi mnobaînoiaiïfëo >b ;iir,rig;2‘)l)'l : *d 

dans lequel a b, — - serait encore un. septième à, p,^ { près de 

a b : d’un autre côté, il résulte des expressions de 0^4-et 
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(n° 3 ) que les termes de Q. qui renfermeront des sinus ou 
cosinus de jabnt , auront en même temps pour facteur la 
puissance j due., ce qui les rendra beaucoup moindres que 
les autres; la plus grande valeur de à’n, résultant de la for¬ 
mule précédente, aura donc lieu dans le cas dey = o ety'—o, 
sans que i soit nul ; et elle sera à très-peu près 

$'n — ^^7 COS. [(a — a) t g' — g)]. 

La valeur correspondante de ja- qui s’en déduira par une 
nouvelle intégration., sera 

*> = c . Sin. [(«' - a) t + 6' — e] ; 

quantité qui ne pourrait être sensible qu autant que * —a 
serait du même ordre de grandeur que le produit HH', ou 
de l’ordre du carré des forces perturbatrices; ce qui ne se 
rencontre pas dans la nature. 

(n) La même analyse s’applique à la partie R de la formule 
( îo). Les développements des quantités P et P r auront la même 
forme que celui de dont ils se déduiront par des diffé¬ 
rentiations relatives à p, rc, et si l’on ne con¬ 
serve dans R que les termes qui ne sont pas à courte pé¬ 
riode, il faudra prendre pour P et P' des termes dépendants 
d’un même multiple de nt. Soit donc 

P=Fcos. {in t-hjabnt — v t —6), 

P* = F'cos. {int+j'abnt — a! t —g'); 

F et F' désignant des coefficients donnés, différents de H et H', 
et les autres constantes étant les mêmes que précédemment. 
On aura 
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i+ja b 


in+f 


~J ab \ 
abn — cl' J 


—-FF '( - i+ j' ab 

\in-\-jabn —a 

COS. [(/—/) a b n t -f- («—«'/ + g — g'], 


en rejetant le terme dépendant de ü n t dont la période serait 
d a peu près un demi-jour. En multipliant par une constante 
y , on aura la partie correspondante Je R ; d’où l’on conclura 
ensuite 


n = T *?„. 1 sin. [(/—/) abnt+{a—a)t + e — ë’], 
oii l’on a fait pour abréger 


[a i + (J+j')a!’} (/—y) + ^ (« — «')+ ! 

Çi+jai- ï) (i + j'a [r/—/> « -I- i («—«')] _ 

Si l’on a /=o et j'=o y les coëffîcients F et F 7 seront in¬ 
dépendants de e ou de/et/*' ; mais si i est nul en même temps, 
cette valeur de S 1 n s’évanouira, et si l’on n’a pas i= o, on 
aura à très-peu près 

i’ n = sin. [O — a') t + ë — ë'] ; 

quantité du second ordre par rapport aux forces perturba¬ 
trices, comme avant l’intégration. La différence «—a'n’étant 
pas de cet ordre de petitesse, la valeur correspondante de 
sera encore insensible. 

Si l’on ay = i et/= t, et qu’on prenne 

p_p/_ dç\ 

df' r —df" 

les coefficients F et F' seront aussi indépendants de /et f. 
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On devra prendre en même temps 

Tl 

^ À ’ 

par conséquent on aura 

^ n F F '(i -b a b) sin. [(g — «') t + € — 6'] 

2AB abnÇi-+-ah — Çi+ab ——^ 

Lorsque i ne sera pas nul, cette quantité sera du'second 
ordre, et son intégrale ou la valeur de £ 2 p.,ne pouvant s’abais¬ 
ser qu’au premier, sera insensible comme précédemment. 
Dans le cas de ^=0, nous aurons 

_ F F' sin. [(« — g.')t-hG —g']~ 

2 AB/2 Ça b — ^ Ça b — — ^ 

et en même temps 

^ 2 FF'cos. [(g—«') t+ g — g'] 

îABoQii— (ab— i) («'—«) 

En supposant le diviseur a' — a de l’ordre des forces pertur¬ 
batrices, on pourrait croire que la petitesse des deux autres 

diviseurs ab — ^ et b —^ suffirait pour rendre sensible cette 

dernière quantité ; mais il n’en est pas ainsi parce que, dans 
le cas de i=o , les coefficients F et F' sont tout-à-fait insen¬ 
sibles, comme on le verra bientôt. 

Cette discussion et celle du n° précédent prouvent donc 
qu’aucun terme de la formule (10) ne peut s’abaisser au pre¬ 
mier ordre par une intégration , ni de deux ordres par deux 
intégrations, et qu’il n’y a non plus aucune circonstance 
numérique qui puisse, malgré cela, rendre sensibles quel- 
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ques termes de En rapprochant cette conséquence de 
celle du n° 8 , il en résulte que la quantité n doit être re¬ 
gardée comme invariable, et la quantité (/. comme propor¬ 
tionnelle au temps ou égale à n t. 

(la) Pour déterminer les variations de fetf, d’où dépen¬ 
dent celles des pôles à la surface de la terre, il est nécessaire 
d’effectuer le développement de Ù que nous avons seulement 
supposé jusqu’à présent. Soit donc L la masse de l’un des 
astres qui agissent sur le sphéroïde terrestre; a, g, y les coor¬ 
données du centre de gravité de L, rapportées respectivement 
aux axes des moments d’inertie A,B, C; .z 0 j,, les trois 
coordonnées relatives à ces mêmes axes, d’un point quel¬ 
conque de la terre; et dm l’élément de sa masse. L’attraction 
étant en raison inverse du carré des distances,, la partie de 
a relative à l’astre L, aura pour expression • 


Ht 


h dii 


[(« - r,T+(f— : 


en étendant l'intégrale à la masse entière de la terre. Chaque 
astre que l’on voudra considérer donnera lieu à un terme 
semblable. Pour avoir égard à sa non-sphéricité, il faudrait 
remplacer L par l’élément de sa masse, et intégrer dans 
toute son étendué; mais il est facile de s’assurer que le terme 
dont le précédent se trouverait augmenté, peut être négligé 
à l’égard de celui-ci, l’un étant à l’autre dans le rapport du 
carré du rayon de l’astre au carré de sa distance à la terre. 
Dans la question qui nous occupe, les astres que l’on con¬ 
sidère sont le soleil à cause de la grandeur de sa masse, 
et la lune à raison de sa proximité de la terre. 

D’après les propriétés du centre de gravité et des axes 

5 
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principaux , on a 

J^x l dm = o^ Jy 1 dmz==yo, Jz,dni^= o; 

Jx t y,dm = o, Jx,z,dm = o, Jy, z,d?n — 0. 

Si donc On appelle £ et r, les distantes deTastfe L et de i’ele- 
ment dm au centre de gravité de la terre, en sorte que nous 
ayons 

ol -h fr‘-i y a — & 1 ', x ? -*-'y? + : « I * = /* x a -, 
et que l’on, djéveloppe suivant lps puissances de #*,/., z,, 
la quantité .comprise sous le signe* J, on aura 

Ù ~ f r;2 ^ m 

-+- ^ Çyjx' dm è J y* dm + y 2 Jz 1 i dm'j-t- D, 

en désignant par D la partie de il qui résulterait des termes 
dp;t/x>iè)dépensions et au-dessus par rapport à x^y, y z,. Les 
couiidoi^ipus dépendront des. angles 6, f, ç, et par 

suite des quantités n, e, c, \yg, l; la distance S en sera in¬ 
dépendante; et comme on n’a besoin que des différences par- 
tieXles de ü relatives à ces quantités, on pourra supprimer 
les,termes de sa,,valeur qui ne contiennent que L Mettant 
dppç <ü l —j a? -H - £ 2 4 \a *pl^ce de y 2 ,, et:observant qu’on a (n° 5) 

fa#/ — z ') dm = C -T- A, dm=^C — B, 

il eiï résultera 


AUTOUR DE SON CENTRE DE GRAVIT^. 3l 

Supposons que L soit la masse du soleil, et désignons par 
m et p sa vitesse angulaire et sa distance moyenne, de sorte 
qu’en négligeant la masse de la terre par rapport à L, on 

ait • ' ' ; - 1 ' : 



Désignons aussi, relativement à la lune, par L', p', 

les quantités analogues à L, p, soit le rapport des 

fractions et ~r dui mesurent lés actions de la lune et du 

soleil, perturbatrices du mouvement de la terre, ou autre¬ 
ment dit, faisons 



En ayant égard en même temps T cés ; deux forces, la valeur 
complète de Q. sera 


-= S 4[(C-A)(¥ + 




+ (C-B)(^£ + =££)]+d. 


(«0 


(i 3 ) La quantité D se composera maintenant de deux 
séries très-convergentes, l’une relative à l’action du soleil et 
l’autre à l’action de la lune. Pour en simplifier le calcul, nous 
considérerons la terre comme un sphéroïde recouvert en 
très-grande partie d une couche fluide en équilibre lorsqu’on 
a seulement égard à la pesanteur et à la force centrifuge; et 
nous ferons usage de la formule connue qui se rapporte aux 
attractions de semblables corps (*). 


(*) Mécanique céleste, tome II, page io3. 


5 . 
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Supposons donc que le rayon r 0 dont l’origine est au 
centre de gravité de la terre, appartienne à un point de sa 
surface, dont v x et u x seront la longitude et la latitude; dé¬ 
signons par r le rayon moyen de la terre, et faisons 

r\ = r'(i + Y,+ Y 3 -h Y 4 + .-t- Y, + etc.); 


Y, étant une fonction rationnelle, entière et du degré i des 
trois quantités cos. sin. w, sin. , sin. w, cos.i;,, qui satis¬ 
fait en outre à cette équation aux différences partielles : 


flLsin 


dYi 


d' Y, 


sin.’a d'v J 


+ i(i 4- i) Y. = o. 


Supposons aussi le centre du soleil situé sur le prolongement 
du rayon de manière que ses trois coordonnées a, 6, y, 
aient pour valeurs 


*.= £sin. u t sin. , g = £sin. w,cos. v,, y= 8 cos.u l . 

Enfin désignons jpar M la masse de la terre, et par 7 le rap¬ 
port de la force centrifuge à la pesanteur sous l’équateur. 
D’après la formule citée, on aura 


- ML MLr- 

n =— + -T- 


[ixteos.-tt, —i) + Y, + ÇY, 


'T Y< - 


■ etc.J , 


pour la partie de fi relative à l’action du soleil. 

Si l’on met dans la première valeur de fi en série du n" 
précédent,à la place de a, g, y, leurs valeurs, et qu’on la com¬ 
pare ensuite à celle-ci, on en concfura 
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(3sin. I w, sin.*^— i) Jx* dm -h ^-^sin.’w.cos. 1 v t — i) j'y'dm 
+ f^cos.’w,— i) j'z, 1 dm=Mr 1 [j^(cos. 2 u, — j ) h- Y,], 
et en outre 

D =—3^— (Y, + y Y 4 -Ketc.)- 

La partie de D relative à l’action de la lune se déduira de 
cette dernière formule, en y changeant L et £ en L' et et 
supposant que les angles u t et v t répondent à la direction 
de son rayon vecteur. 

Si la terre était elliptique, la quantité Y, subsisterait seule 
dans l’expression de son rayon variable, et toutes les sui¬ 
vantes Y 3 , Y 4 , etc., seraient milles. L’ensemble des nom¬ 
breuses mesures des degrés et du pendule que l’on a faites, 
s’écartant fort peu de cette hypothèse, il en faut conclure 
que Y 3 , Y 4 ,etc., sont,sinon égales à zéro, du moins très-pe¬ 
tites par rapporta Y,. D’ailleurs elles sont multipliées dansa, 

une fois, deux fois, etc., de plus que Y, par j ou ^ , c’est-à- 

dire par la parallaxe du soleil ou de la lune; il en résulte 
donc, pour cette double raison, que la partie D de la for¬ 
mule (il) sera très-petite relativement à l’autre partie, dé¬ 
pendante de C — A et C — B, qui formera à très-peu près la 
valeur de Cl. Néanmoins nous conserverons, pour plus d’exac¬ 
titude, la partie D réduite à son premier terme. 

En regardant la terre comme un solide de révolution, la 
quantité Y 3 que ce ternie contient, aura pour expression : 

Y,=SGcos.’«,-4cos .».)=5 (£t-&); 
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i étant une constante dépendante de la différence d’apla¬ 
tissement des deux hémisphères, qu’on devra supposer peu 
considérable par rapport à l’aplatissement moyen, et, par 
exemple, au-dessous d’un dix-millième. Faisons en outre 

ü 5 <7 ’ 


n étant une quantité qui diffère peu de l’unité d’après l’ob¬ 
servation. ; Si l’on forme les parties de D relatives aux actions 
du soleil et de la lune; que l’on élimine les masses L et L' 
comme précédemment, et ia masse M au moyen de cette der¬ 
nière équation , on aura 



4 Cm 




P 3 Y 3 


5 r' p 3 f wr' p /a y' 3 

ï¥ + F 


5 wr'p' 3 y' 

~~W ~ 5 


)■ (r3) 


Cette valeur de D suffira, en général, à l’objet que nous 
nous proposons; elle ne contient, comme on voit, que des 
puissances impaires de y et y' ; elle renfermerait des puis¬ 
sances semblables des autres coordônnées a, a', 6,6', si 
l’on avait égard à 1 a non-symétrie de la terre autour de son 
axe : nous donnerons plus Las un exemple du peu d’influence 
sur le mouvement de la terre, des termes de D qui résulte¬ 
raient de cette considération. 

(i 4 ) La valeur la plus générale de Y 2 est de la forme : 

Y, = ■/}(-£ — cos. 2 -h7i / sin. 2 n l cos.2v 1 + ri"sin. 2 u x sin.2v t 

+ /'sin. 2M, cos. v, -+- 7 ) iy sin. 2 u t sin. v t ; 

y",?T, étant cinq coefficients constants. En la sub¬ 
stituant dans l’équation (12), on en conclura en premier lieu 

y," = o, y) , " = o, v) ,r = o ; 
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ce qui provient de ce que les droites d’où les angles u, et v, 
sont comptés, sont des axes principaux. Si L’on exprime en¬ 
suite les intégrales que cette équation renferme, au moyen 
des moments d’inertie A, B, C (n° 5 ), on trouvera quelle peut 
s’écrire ainsi : 

|(2C—A — B){| — cos.’m,) -h 7(A—B) sin/w, cos. 2v t 

=Mr h (n — 7x)(t — cos. 1 ) + Mr'* V sin. 1 w, cos. 2 v t . 


Elle se décomposera donc en deux autres, savoir : 

A — B=| M/’n', 

2C—A — B = | Mr' , (n — {%). 

L’aplatissement d’un méridien quelconque, résultant de 
la valeur de Y,, sera exprimé par n 4-r/co$. nv t . D’après les 
observations du pendule, il est sensiblement le même pour 
tous les méridiens ; il faut donc que V soit très-petit par rap¬ 
port à », et aussi par rapport à *—7^, parce que n’est pas 
le tiers de » ; par conséquent la différence B—A est très-pe¬ 
tite relativement à 2 C — A—B, et aussi à l’égard de chacune 
des différences C—A et C — B. 

La dernière équation, mise sous la forme 


2 C —A —B 
C 


_ u ( _ _, 




(>4). 


au moyen de la valeur de C du n° précédent, nous sera utile 
dans la suite de ce Mémoire. 


(i 5 ) Cela posé, si nous désignons par 2, les trois 
coordonnées rectangulaires du centre du soleiL, relatives à 
des axes menés par le centre de la terre suivant des direc¬ 
tions fixes; que nous prenions pour celui des x y la droite 
d’où l’on compte l’angle et pour celui de z la droite per- 
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pencLculaire au plan de cet angle (n° i), nous aurons, d’après 
les formules de la transformation des coordonnées, 
a = x (cos. G sin. 4 sin. 9 -+- cos. 4 cos. 9) 

4 -/(cos.G cos. 4 sin.9 — sin.^cos. 9)—■ 5;sin. G sin. 9, 
ê =x (cos. G sin. 4 cos. 9 — cos. 4 sin. 9) 

+ /(cos. G cos.4 cos. 9 h- sin. 4 sin. 9)—zsin.Gcos. 9, 
y = #sin. Gsin. 4 + /sin. G cos. 4 -t-zcos. G. 

On en déduira a',ê',y', en y remplaçant a?,/,z, par les 
coordonnées analogues du centre de la lune que nous repré¬ 
senterons par #',/', z'. 

Nous supprimerons, comme plus haut, dans les différences 
partielles de 12, les termes dépendants de l’angle nt-±- 1 . Or, 
en considérant les valeurs de G, 4 et 9 données par les équa¬ 
tions (4), il est facile de voir que les parties de a 2 , 6% y 3 , y, 
et de semblables quantités relatives à la lune, qui dépen¬ 
dront de la première puissance de e, ne contiendront que 
des puissances impaires de nt+l ; il en sera de même à 
l’égard de la partie de 12, linéaire par rapport kf et/", et 
résultant des formules (11) et (1 3 ) ; donc aussi les valeurs 

correspondantes de ne renfermeront que des ter¬ 

mes à courte période, et devront être entièrement suppri¬ 
mées. Ainsi en négligeant les termes qui ont pour facteur /ou 
/', les équations (9) se réduiront à 

df= o df' = o; 

et dans cette première approximation, les quantités/'et f 
seront constantes abstraction faite des termes insensibles 
dont la période est un jour ou un sous-multiple du jour. 

(16) Si l’on tenait compte de la non-symétrie de la terre 
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autour de son axe, cette circonstance introduirait dans D 
des puissances impaires de a,a',6,6', qui donneraient lieu 
a des termes indépendants de n t + / parmi ceux du premier 
ordre par rapport a e. Il est bon d’en apprécier l’influence, 
et de montrer quelle serait tout-à-fait insensible. 

Pour cela, soit 

psin. w, cos. 

un des ternies de Y 3 dus a cette cause; p étant un coefficient 
constant, tres-petit par rapport à l’aplatissement de la terre, 
et, p?.r exemple, au-dessous d’un dix-millième. La partie 
correspondante de D ou de ü, sera 

en faisant, pour abréger, 

5 g r' 5 w g r' 

2p 2 p' S 5 

et d après les quantités qui entrent dans s et s', on peut les 
évaluer avec une approximation suffisante, à 

^ = 0,0001, s' = o,i. 

En mettant pour a et a.’ leurs valeurs, après y avoir sub¬ 
stitué celles de 0, <|> et 9, on trouvera que la partie de n 
linéaire en/’et f\ et indépendante de l’angle /, se com¬ 

posera de termes de cette forme : 

fc^WN/+ N'/'jcos.^ènî-vf-c,) 

+ s (Qf+ Q/ )cos. {abnt —v' t —<?')]; 

les coefficients constants N, N', Q , étant des fractions peu 

6 
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considérables à cause des facteurs a ou b qu’ils renfermeront; 
v,v', c l et c', désignant aussi des constantes données, et vt 
provenant du mouvement du soleil, et v ' t de celui de la lune. 
En faisant en dehors des cosinus A = B = C et abn = m, 
ce qui est suffisamment exact, la première équation (g) don¬ 
nera 

df— — [.y N' cos.(ab n t —v£—c,) 

-h s Q'cos .{abnt — V t —c')]^/?* ndt r 

et en intégrant 

f— p, n N'sin. ( abnt — vt — c t ) 

-f- ■ m ,— Q' sin .(abnt — V t — c)\ 

La valeur de f\ donnée par la seconde équation (9), se dé¬ 
duira de celle de f , en changeant le signe, et y mettant N 
et Q à la place de N' et Q'. Il suffira d’examiner la valeur 
de/ 

Le diviseur v — abn aura la plus petite valeur quand vt 
sera le moyen mouvement du soleil; il sera alors égal à en¬ 
viron un septième de m, et l’on aura 


Si Vt est le moyen mouvement de la lune, ou treize fois 
celui du soleil, il en résultera 

s'm 1 

—-7- <-. 

v —abn 120 

Cette dernière quantité sera à peu près douze fois plus grande, 
lorsque V t sera le moyen mouvement du périgée ou du nœud 
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de la lune; mais dans ce cas, le coefficient Q' par lequel elle 
est multipliée aura pour facteur l’excentricité ou l’inclinaison 
de l’orbite lunaire sur l’écliptique. Tout cela joint à la peti 1 - 

tesse du facteur p, nous fait voir que le rapport résultant 

de la formule précédente , ne sera qu’une fraction insensible 
qui ne s’élèvera pas à un millionième. Il en sera de même à 

1 egard de — ; d’ou nous concluons que la non symétrie de 

la terre autour de son axe ne fera pas varier sensiblement 
les quantités/et/; ce qu’il s’agissait de démontrer. 

D’après la valeur précédente de a, les coëfficients que 
nous avons représentés par F et F' à la fin du n° 11, auront 
pour expressions : 

F = ^> icN r== ^> lcN 
n ' n 1 

ce qui justifie ce que nous avons dit sur le degré de petitesse 
de ces quantités et du terme de 8 1 ^ qui dépend de leur produit. 

(17) Si nous voulons connaître les inégalités séculaires 
qui peuvent affecter les valeurs de/et/', il faudra conserver 
dans les seconds membres des équations (9) les termes qui 
renferment les premières puissances de ces quantités, et 
dont les coëfficients sont indépendants, non-seulement de 
nt + /, mais encore de l’angle a b (nt •+- c), des moyens mou¬ 
vements du soleil et de la lune, et de ceux du périgée et du 
nœud de l’orbite lunaire, afin qu’on puisse les considérer 
comme constants. Les termes qui renferment n indépendam¬ 
ment de nt dans les valeurs de 6, ^ et ç ayant pour facteur 

e, il en sera de même à l’égard de par conséquent les 

6. 
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seconds termes des formules (9) seront du deuxième ordre 
par rapport à fetf\ et, comme tels, devront être négligés. 
On ne conservera donc que les premiers qui de'pendront de 
la partie de H dépendante du carre' de e. 

Pour l’obtenir, prenons le plan de l’écliptique à une époque 
déterminée pour celui des j, sur lequel on compte l’angle 
+ ; ce qui permettra de supprimer l’ordonnée z, dans les va¬ 
leurs de a, 6, y. Si nous négligeons, en outre, l’excentricité 
de l’orbite solaire, et que nous représentions par v la lon¬ 
gitude du soleil, comptée de la même origine que mais 
en sens contraire de cet angle, nous aurons 

£=:p, ^ = pcos.'i;, j=rsin.o>, z= o; • 

d’où il résultera (n° i5) 

a=p[sin.( / u 4- <J/)cos. 9 sin. 9 -h cos. (v 4- <|/)cos. 9,], 

6 = p[sin.('y 4- <]/)cos. 9 cos. 9—cos.(i> -h ^)sin.(p], 
y === p si n. ( 4 - <j>) sin. G ; 

formules dans lesquelles v 4- <|> sera la longitude du soleil, 
comptée du nœud descendant de l’équateur. En rejetant les 
termes dépendants de cet angle, on aura simplement 

a 2 =7 p 2 (cos . 2 6 cos.*9 4- sin . 2 9) ,- 
ê 2 =7 p 2 (cos . 2 9 sin. 2 9 4- cos. 2 9); 

et Tes puissances impaires de y qui entrent dans la formule 
(i 3 ) seront égales à zéro. 

Je substitue dans a 2 pour 9 et 9 leurs valeurs; je néglige 
les puissances de e supérieures au carré, ainsi que les termes 
dépendants de l’angle ab(nt-hc), ce qui fera disparaître les 
premières puissances de k et À'; de cette manière, il vient 



AUTOUR DE SON CENTRE DE GRAVIT^. 


43 


* ,== 7 p* [cos. 2 x cos. * (11 1 1) 4 -sin 2 («£ 4- /)] 

p a e a r 

— ~Tn 7 K iCOS -’* + t^'cos. sx h- [A'(i 4- cos. 2 x)—2Â*X'cos. 2 x]sin. 2(« t-\- /) 
4- (k 1 * cos.*X — h COS. 2 X— k 1 COS. 2X) cos. 2 (/l £ 4- /)j ; 
et l’on aura en même temps 

£= g^[ B ( G — A ) 4 - A ( c — B ) — C ( B —A)cos.2(rcf 4 - /)], 
h. =-g^-sin. 2 (wf+/)» 

^ ,= = — A) 4 - A(C — B) 4- C(B— A)cos .^(nt 4 - /)], 

^ ,=:: ^-C^[B(C— A) + A (C—B) — C(B—A) cos. 2 (nt 4 - /)], 
X'^'=^*[B(C—A) 4 - A(C — B)]sin. 2 (nt 4- /); 

au moyen de quoi, Ton trouve 

2 A (G — B) p a c*sin.*X 

“ 8C^ » 

en supprimant la partie indépendante de e, qui disparaîtrait 
par les différentiations relatives à fe t/', et négligeant les 
termes dépendants de ni 4- L On trouvera de même 

_B(C — A) p’e a sin. a X 

“ 8G 

Si l’on néglige l’excentricité et l’inclinaison de l’orbite 
lunaire, les valeurs de a ' 2 et ê ' 2 se déduiront des précédentes 
en y changeant p en p'; l’ordonnée y' sera nulle ainsi que y, 
et par suite la valeur de D donnée par la formule (i 3 ). En 
vertu de l’équation (11), on aura donc 
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fl = 3OT ’ ( , 1 6^ Sin ' n ( A + B ) (C- A )(C-B)(/‘ +/'•). 
Par conséquent, les équations (9) deviendront 
df— — A mfdty df l = t^mfdt, 
en faisant pour abréger, 


3 m 
“ 8 n 


(A + B ) ( C-A )(c -B) (l+M)8in , x=A; 


et si nous désignons par E et T les deux constantes arbi¬ 
traires , nous aurons en intégrant : 

f— n E cos. (Amf + r), f'=znTL sin. (A m t r ). 

D’après les données qu’on a sur les quantités qui entrent 
dans A, la période de ces inégalités séculaires de fe t f sur¬ 
passerait cent millions d’années : si elles étaient sensibles, ce 
seraient les plus lentes de toutes celles qui existent dans le 
système planétaire. 

(18) Les valeurs correspondantes de données par 

les équations (4) seront 

pz=na^QOS.{abnt —A mt —r), 
q =znbYLs\n;iabnt — A mt — r), 

r — n\ i — -- - -t--——cos. liabnt — Amt —r)J* 

A cause que abn est un peu plus grand que m , la période 
des valeurs de p et q sera d’un peu moins d’une année. En 
désignant par A la distance angulaire de l’axe de rotation à 
celui du plus grand moment d’inertie C, on aura 

sin. A — v/ — = a E , 

^ + ^ + 
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en négligeant le cube de E, et observant que l»es quantités 
cl et b sont.a très-peu près égales. Le pôle de rotation décri¬ 
rait donc chaque année à la surface de la terre un cercle d’un 
rayon égal à «E, celui de la terre étant pris pour unité; 
d où il résulterait, dans les latitudes géographiques, une va¬ 
riation annuelle double de aE; or, l’observation n’ayant fait 
reconnaître aucune variation sensible dans la latitude d’un 
même lieu , on en peut conclure que l’angle aE est au-des¬ 
sous d’une seconde, et en prenant^ pour la valeur de a, 
d s’ensuit que E est au-dessous d’un dix-millième. Mais il y 
a heu de penser que cette limite est èncore trop élevée, et 
que le coefficient a E est du même ordre de petitesse que les 
forces perturbatrices, c’est-à-dire du même ordre que le pro¬ 
duit du rapport ^5-et de l’aplatissement de la terre. 

Les termes du développement de ü qui dépendront de l’angle 
abnt, ayant pour facteur cette quantité aE, on vérifie mam- 
tenant ce que nous avons dit précédemment (n° 8), savoir, 
qu’ils ne peuvent avoir aucune influence sensible sur l’inté¬ 
grale seconde de la formule (8), lors même qu’ils acquer¬ 
ront par une double intégration le diviseur (ab — , ce 

qui ne les empêchera pas de rester au-dessous de l’ordre 
de grandeur des forces perturbatrices. 

La vitesse de rotation 8 déduite de ces valeurs de p, q,r 
(n° 1), aura pour expression : 


*[■ 


(A -f- B) (G—A) (C — B) E ! 
4 ABC 


EV B—A , , 

V 1 + J^x^s.2(abnt—Kmt— T 


Donc, d’après ce qui a été démontré à l’égard de * et de Fin 
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tégrale j^ndt, les variations de 8 seront insensibles, et l’in¬ 
tégrale j'tidt sensiblement proportionnelle au temps. 

Nous tirerons de tout ce qui précède les conclusions sui¬ 
vantes : 

i° L’action du soleil et de la lune sur le sphéroïde 
terrestre ne produit aucun déplacement des pôles de rota¬ 
tion à la surface de la terre, et l’axe de rotation coïncide 
constamment avec celui du plus petit moment d’inertie; 
d’où il résulte que la longitude et la latitude de chaque lieu de 
la terre sont invariables, ce qui n’aurait pas lieu si l’axe de 
rotation ne rencontrait pas toujours sa surface dans les mêmes 
points. 

2° La vitesse de rotation est aussi sensiblement invariable. 

3 ° L’intégrale de cette vitesse multipliée par l’élément du 
temps, ou l’angle décrit par chaque point de la terre en 
vertu de sa rotation, croît proportionnellement au temps 
employé à le décrire, abstraction faite de petites inégalités 
qui seront toujours insensibles. 

Les deux premières propositions étaient déjà démontrées 
dans le V e livre de la Mécanique céleste et dans mon Mé¬ 
moire sur la Rotation de la Terre . C’est dans ce Mémoire que 
la troisième a été démontrée pour la première fois, par une 
analyse différente de la précédente. Elle n’est pas une con¬ 
séquence nécessaire de la seconde ; car pour parvenir à celle- 
ci, il suffit, en général, de considérer les termes du premier 
ordre par rapport aux forces perturbatrices, et pour dé¬ 
montrer l’autre, il était nécessaire d’avoir égard , comme nous 
l’avons fait, aux termes du secondordre. 
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L’intervalle de temps pendant lequel l’intégrale jttdt 

augmente de quatre angles droits, est la durée du jour si¬ 
déral; cette durée est donc constante, et les jours sidéraux 
successifs peuvent servir de mesure au temps, ainsi que 
cela se pratique effectivement dans les observations astrono¬ 
miques. 

(19) Puisque les quantités ae et be sont regardées comme 
insensibles, il en sera de même à l’égard de k,k\'h,h\ et 
les valeurs de G, ÿ et 9 du n° 3 , se réduiront à 

6=A, ty=g, <pt=znt+l. 

La fonction û. ne contiendra plus n indépendamment de 
nt + l; par conséquent la différence partielle ~ qui entre 

dans la valeur de dl sera nulle; et en vertu des deux der¬ 
nières équations (7), on aura 

dl=cos.\dg; 

ce qui fera connaître l au moyen de g ou de Mais il suffira 
de déterminer les angles G et ÿ en fonctions de t, pour fixer 
à chaque instant la position de l’équateur terrestre et de 
l’axe de rotation qui lui est perpendiculaire. 

Si l’on néglige les quantités du second ordre par rapport 
aux forces perturbatrices, la différentiation relative à / ne 
laissera subsister dans ^ que des termes dépendants de 

qui n’augmenteront pas par l’intégration; en les sup¬ 
primant donc, désignant par V la partie de il indépendante 
de nt+ et mettant G et ^ à la place de 1 et g, les quatrième 
et cinquième équations (7) deviendront 


7 
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7 dW dt 

dty C rasin.ô 1 
dV dt 
d§ C«sin.0 

Ainsi les déplacements de l’équateur par rapport à des plans 
fixes, se détermineront au moyen de ces deux équations très- 
simples auxquelles j’étais déjà parvenu d’une autre manière. 

On peut remarquer que les coordonnées 6 et 6', a et a' qui 
entrent dansft, se déduisent l’une de l’autre en augmentant <p 
d’un angle droit (n° i 5 ); d'où il résulte que quand on rejet¬ 
tera, pour former V, les termes dépendants de l’angle 9 et 
de ses multiples, on aura ë 2 = a 2 , é' 2 = a' 2 , et, par conséquent 
(n° 12), 

V = ^(aC—A-B)(^+!^^)+D; (i 6) 

la partie D étant toujours donnée par l’équation (i 5 ). 

§ III. 

Formules de la précession et de la nutation. 

(20) Prenons, comme précédemment, le plan de l’éclip¬ 
tique à une époque déterminée pour celui de l’angle , qui sera 
aussi le plan des coordonnées x et ÿ et/; et rappelons- 
nous que ^ sera la longitude du nœud descendant de l’equa- 
teur, comptée à partir d’une ligne fixe, en sens contraire de 
la rotation de la terre (n° 2). Désignons par v la longitude 
de la lune, par V celle du nœud ascendant de son orbite, et 
par t é celle de son périgée; ces trois angles étant comptés à 



AUTOUR DE SON CENTRE DE GRAVITE. 49 

partir de la même ligne fixe que <|>, mais en sens contraire, 
de sorte que?/4-^,^4-<j>,7r' 4-<|/ soient les longitudes comp¬ 
tées, à l’ordinaire, dans le sens de la rotation et à partir 
du nœud descendant de l’équateur. Représentons par l’in¬ 
clinaison de l’orbite lunaire sur le plan fixe, et par e' son 
excentricité. Faisons, pour un moment, 

[ i 4 - tang. a )/sin. , ('i/ — /')] ’=U; 

nous aurons, sans rien négliger, 

x f =Ucos .v , ^' = U&'sin. V, z'=U£'tang.>'sin. {y —/'), 
et par conséquent 

a '=U [cos. Gsin.^' -f- <J;) sin.9 -h cos.^i/ 4- ^)cos. 9 

— sin. 0 tang.Vsin. (V—/')sin.9j. 

Si l’on supprime les termes qui contiendraient sin.^ou 
cos.29, on en déduira 

a ,, = jU , ^ ,ï j I 4- COS.’e 4- sin/Ocos.^^' 4 - «!/) 

4- sin.’G tang.)/[i —cos^o/— /')] 

— 2sin.0cos.0tang.x'[cos.(<|> 4- /')—cos.^?/4- <|>— /')]J. 

Par les formules du mouvement elliptique, on a 

— = 1 +-é* — e cos.fm't 4 - t — r:')—~e' i cos.Q.{m't-\-i '—*'), 
P 

v'=m't+ e '4-2c'sin.(7n'^4- e '—*')+{ e"sin.ï(m t-\-i —x') 
— f>' J sin.2(ttï'£ 4-e’ — /'), 

en désignant par m' 1 4-e' la longitude moyenne de la lune, 
et négligeant les termes du troisième ordre par rapport à 
et é . Au même degré d’approximation, et d’après ces valeurs 

7 - 
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on trouvera 


a' a p' 3 , / 3e ,a .n 

-rX 2 J(r + cos. 2 ô)4-^x' 2 sin. , e 

—7Vsin.0cos.0cos.(«(» + 1 ') + ^'’sin.’Oeos.a(<f -h/') 

+ 7sin. a 0cos. 2(771 1 + e' + : 


les termes périodiques qui dépendent de m't+ e' augmen¬ 
tant beaucoup moins que les autres par l’intégration, on a 
négligé l' et e dans leurs coefficients, ce qui les a réduits au 
dernier terme de cette formule. Si l’on avait égard aux per¬ 
turbations de la lune dans les valeurs de X et v\ il en résul¬ 
terait d’autres termes périodiques dépendants des moyens 
mouvements de la lune et du soleil, qu’on peut négliger rela¬ 
tivement a ce dernier terme. Celui qui proviendrait de la 
'variation et qui dépendrait du double moyen mouvement 
du soleil, croîtrait a la vérité dans un plus grand rapport par 
l’intégration; mais d’après la grandeur de cette inégalité de 
la lune, il serait environ un centième du terme semblable, 
dû à l’action directe du soleil sur le sphéroïde terrestre, et 
l’on peut en conséquence en faire abstraction. Les inégalités 
principales dont \ et V sont affectés, et qui ont la même pé¬ 
riode que la variation, donneraient lieu à des termes encore 
plus petits. 

Relativement au soleil, désignons par mt+& sa longitude 
moyenne, par w celle de son périgée, par e l’excentricité de 
l’orbite, par 1 son inclinaison sur le plan fixe des x, j, et 
par l la longitude du nœud ascendant sur ce même plan; la 
formule précédente s’appliquera immédiatement à cet astre, 
en y supprimant les accents ; de plus l’inclinaison 1 de l’éclip¬ 
tique mobile sur 1 écliptique fixe étant une très-petite quan- 
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tité, on pourra négliger son carré; et alors on aura simple* 
ment 

+ cos.* 6 )—-^sin. Ôcos. 6 cos. (<(/4-/) 

-h ~ sin. 2 ô cos. 2 (m 1 4- e 4 - +). 
Au moyen de ces résultats, l’équation (16) deviendra 
V=}m*( 2 C-A-B)|(i + u + COS- . 

— [xcos. 4 - /) 4 - uV cos. (<|> -h /')]sin. 26 
-h 4- 7 cos. 2 (v|/ 4- r)]sill . a 6 4-[cos. 2(mt 4- e 4-^) 
4 - w cos. 2 {ni 1 4- i 4 - <[>)] sin . 2 61 , 

en faisant abstraction de la quantité D. 

Cette formule ne contient, comme on voit, aucun terme 
dépendant du périgée du soleil ou de celui de la lune. Il n’en 
serait pas de même si l’on tenait compte de la partie D, donnée 
par l’équation (i 3 ); les puissances impaires de y et y' y intro¬ 
duiraient des termes qui dépendraient de 4- tz ou de ^ 4- 
mais en considérant qu’ils auraient pour facteurs l’excentricité 
et la parallaxe du soleil ou de la lune,et qu’en outre le coefficient 
C£ de la formule (i 3 )est beaucoup moindre que le coefficient 
2 C—A—Bdela précédente, on voitqu’ilspeuventêtrenégli- 
ge's par rapport à ceux que nous avons conservés. L’expression 
de D renfermerait aussi des termes dépendants de l’angle 3^4- 
tt' 4- 2 qui croîtraient à pe-u près vingt fois plus par l’intégra¬ 
tion que ceux qui repondent à l’angle ^ 4- r 7 ; mais comme ils 
auraient pour facteur e'\'* au lieu de e', ils seraient encore 
moindres que ceux-ci, et négligeables à plus forte raison. 

(21) L’inclinaison moyenne de lorbite lunaire sur Téclip- 
tique mobile, est, comme on sait, sensiblement invariable; 
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pour cette raison, il est bon de l’introduire dans la formule 
precedente à la place de l’inclinaison sur 1 écliptique fixe. 

Or en la désignant parc, et négligeant les quantités du troi¬ 
sième ordre par rapport à X, a' et c, on aura 

}/sin./ f = csin./' -hXsin l, 

\ 1 cos. V = c cos. V h- \ cos. I. 

En effet, prenons dans le plan de l’orbite lunaire, un point 
dont la distance au soleil soit l’unité; et représentons par x 
sa longitude sur l’écliptique fixe. Ses distances à ce plan et 
à l’écliptique mobile seront exprimées par ^'sin.(a? — V ) et 
csin. (x — l '), au degré d’approximation où nous nous arrê¬ 
tons; de plus la partie de la première distance, interceptée 
entre les deux écliptiques , aura pour valeur xsin.(.æ — /); et 
enfin, cette partie pourra être considérée comme l’excès 
de la première distance sur la seconde. On aura donc 

*sin. (x —/)='x'sin. (x — I) — csin.(a? — /'), 

quel que soit x; et en y faisant successivement x=o et x 
égal à un angle droit, il en résultera les deux équations pré¬ 
cédentes. 

Par leur moyen, et en négligeant le produit \c, la valeur 
de V deviendra définitivement 

V = 4to-( 2 C— A-B) |(i +<* + 3 -f + ^f— ^)(i + cos.'6) 
— [( t + &>) X cos. (+ + /) + « C cos. -h /' )] sin. 2 e 
4-û)C a [i H- 7cos. 2(+ - 4 - 1 ’)] sin. a 6 
-4 -[cos.2(/?zé + e + <]>) -t-tocos. -t- <|;)]sin. a () J ; 

et en la substituant dans les équations (i5), nous aurons 
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db= - n - ~ j[(i 4 ~ tù)xsin. (<|/ 4 - /) -f wcsin. 4 -/')] cos. G 

— [fuc’sin.a^ 4 - /')+ sin. 2 (m£ 4 - é 4 - <j>) 4- «sin. 2 (ra'f 4 - e' 4 -<J/j]sin.G [ dt, 


d 4 


_3/ra* (aC— A—B) 

4«C 


!( 


I 4- a) 4- - 


2 2 J 


cos. 0 


4- [( I 4- COS. (<|> 4- l'j 4* <àC COS. (<J> 4- V )] - '-n 
— [7 (*> C* COS. 2 4- COS. 2, [pit-t -1 4~<J/) 4~ 10 COS. û (jll 1 1-\- e'4-+)]cOS. 0 j dt. 

(aa) Pour intégrer ces expressions, supposons que le temps 
t soit compté de l’époque où l’écliptique mobile coïncidait 
avec le plan de l’angle^, et donnons pour origine à cet 
angle et aux autres longitudes, l’équinoxe du printemps de 
la même époque. Nous aurons à la fois t=o ,x=o, <f=o. 
Désignons à cet instant par e, et A, les valeurs de e et de 6 * 
observons que les variations de e sont insensibles, aussi bien 
que celles dec; négligeons dans une première approximation 
celles de G et de e, la quantité X, et les termes périodiques 
qui dépendent des angles /' ,mt 4 - e, m t 4 - e'; faisons enfin, 
pour abréger, 


y 3 w 2 (2C—A—B)^ . .. 3e* 3(ùe'* 3<oc*\ 

4 ^G { !+*> + — 4-— -—J cos. h. 

Nous aurons 



d G = o, d^z= o, 6 = 4 , <J, = Çr. 

Les inégalités séculaires de e,\sin.l et Xco s.l, n e sont 
pas assez bien connues pour qu’on puisse en faire usage dans 
les valeurs de dü et dty; mais à cause de leur extrême len¬ 
teur, on peut les supposer développées suivant les puissan¬ 
ces de t, et s’arrêter à la première. Soit donc 
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e = e 1 ~b/?, *sin./=£*?, xcos .l=gt; 

f,g->g\ étant des coëfficients constants qui dépendent des 
masses des planètes. Substituons ces valeurs et les précé¬ 
dentes de 6 et ^, dans les expressions qu’il s’agit d’intégrer; 
négligeons le carré de t et les inégalités périodiques ; il en 
résultera 


rf(= i,’( a C-A- B )(, + .) eo i i . 

4«C ° 

— - COS, ^ [3e,/’-t-2(i +u>)g' cot.2 K\tdt; 


et en ayant égard à ce que représente Ç, et intégrant, on en 
conclura 

6 = h - 4 - 7 Çg t 2 , 

( a(^ 8 +âj +g' COt-z/i^Çt'. 

Quant aux termes périodiques dont nous avons fait ab¬ 
straction, il suffira de les intégrer en regardant 6 comme 
constant et <J> comme nul; on pourra aussi n’avoir égard qu’à 
la partie constante du moyen mouvement du nœud de la 
lune, de sorte que si l’on fait 



V — Çj t -+- 6 , 

g' et 6' seront des quantités constantes; et vu la petitesse des 
termes qui renferment mt-\- e et m't-t- e, nous y remplace¬ 
rons ces angles par v et après l’intégration. En appelant 
alors 0 et Y les parties périodiques qu’il faudra ajouter aux 
formules (18), nous aurons 
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®--(7^ cos -( ef+ê 0+^frcos. a(ê(+ r) 

£ tan g. 4 / cos. 2 v cocos. 2 2/\ 

^~a(i + to)V m m! J 1 

T = . sin (g t + e>)__?gL g sin. a(6 f + O 


2 (i + w ) 


/*sin. 2 y «sin. 2 V'\ 
\ m m' ) 



On devrait ajouter à ces valeurs des constantes qui les ren¬ 
dissent nulles quand t = o ; mais on peut s’en dispenser en 
supposant que h représente maintenant l’inclinaison ô qui 
a lieu à cette époque, diminuée de la valeur correspondante 
de 0, et que l’origine de l’angle <|j est lequinoxe de la même 
époque, corrigée de même de la valeur de 

(s 3 ) Les formules (18) déterminent à chaque instant la 
position moyenne de l’équateur par rapport à l’écliptique 
fixe. Elles ne different de celles de la Mécanique céleste (*) 
qu’en ce que nous avons conservé dans la seconde le très- 
petit terme dépendant de la variation de l’excentricité de 
l’orbite solaire, et développé les inégalités séculaires suivant 
les puissances du temps. On s’assurera de leur coïncidence en 
développant ainsi les formules de M. Laplace, jusqu’à la 
seconde puissance inclusivement, et observant que, d’après 
ses notations comparées aux nôtres, on a 

2csin.6=o, 2ccos.ë=o, 2 c/cos. G=g, 2c/sin.é=— g. 

Soit actuellement 0' l’inclinaison de l’équateur sur l’éclip¬ 
tique mobile. Projetons l’équinoxe du printemps sur l’éclip- 


(*) Tome II, page 3i8. 


8 
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tique fixe, et désignons par <]/ la distance angulaire de cette 
projection à la ligne fixe d’où l’on compte l’angle , en sorte 
que <|> — <J' soit sa distance au nœud descendant de l’équateur 
sur l’écliptique fixe. Pour déterminer 8' et concevons une 
sphère dont le centre soit celui de la terre, et considérons 
le petit triangle sphérique formé par les trois arcs compris 
entre les intersections de l’équateur et des deux écliptiques. 
Le côté compris sur l’écliptique fixe sera <], + /; et comme la 
longitude l appartient au nœud ascendant de l’écliptique 
mobile sur ce plan, il en résulte que l’angle opposé sera le 
supplément de 8 / . Les deux autres seront 8 et}.; on aura donc 

cos. 8 ' = cos. 8 cos. \ — cos. (<J/ -+- /) sin, 8 sin. x. 

Appelons l le côté compris sur l’équateur ; abaissons de l’équi¬ 
noxe du printemps, un arc perpendiculaire au côté op¬ 
posé; le segment adjacent à l’équateur sera ^ — <|/, et dans le 
triangle rectangle dont il fait partie, on aura 

tang. (<|>— <]/')=tang.£cos. 8. 

Enfin dans le premier triangle, nous aurons aussi 

. . sin. (ib-f- /)sin. "k 

sin. $ =- VY . - 

sin. 8 

Ces trois équations sont rigoureuses; mais en les résolvant, 
nous négligerons les puissances dex supérieures au carré. La 
première donne alors 

8 =8 -t- \ cos. (tj; + /) + 7 V sin.*(^ +cot. 8 ; 
on tiré de la troisième 

^_’ Xsin.^-f- IJ X*"sin.a (i|> + /)cc>s.8 

sin. ô asin.’Q ’ 
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et ensuite de la seconde 

—Xsin.(<|> -h/)cot. 0 + ^Vsin.(^ -t-/) cot. 3 ô. 

On a du conserver le carré de X dans ces valeurs de ô' et 
'i' , puisque nous poussons l’approximation jusqu’aux termes 
dépendants du carre du temps. Par la même raison, les va¬ 
leurs de Xsin. / et Xco s. I relatives à l’écliptique moyenne 
que nous y substituerons, seront de la forme : 

xsin .l=gt+kt', xcos./=^' 

g et g étant les mêmes que précédemment, et k et k' d’au¬ 
tres coefficients constants. Au moyen des équations (18), on 
aura de cette manière 

b' = h +g't—(igt—±g>cot./i — k')r, 

+ — (ï—£cotÀ)—+ -k-j.—gg-cot.’h + *cot ;.h) t\ 

( 24 ) Ces valeurs de Ô' et f se rapportent aux positions 
moyennes de 1 equateur et de l’écliptique; pour avoir les 
valeurs de ces angles qui répondent à la position vraie de 
l’equateur, il y faudra ajouter celles de 0 et W qui sont don¬ 
nées par les équations (19); et pour quelles répondent aussi 
à la position vraie de l’écliptique, on y ajoutera encore d’au¬ 
tres termes provenant des inégalités périodiques du soleil 
en latitude. Mais dans la pratique de l’astronomie, on réduit 
les positions observées de cet astre à l’écliptique moyenne ; 
ce qui dispense de faire subir cette seconde correction aux 
valeurs de ô' et y. Toutefois, il ne sera pas inutile de former 
les expressions complètes de ces angles, tels qu’ils résulte¬ 
raient de l’observation immédiate. 


(20) 


8. 
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Soit donc z la latitude du soleil; de sorte qu’on ait 
z = Xsin. (y — l ), 

en négligeant toujours le cube de X. Dans la différentielle 
première de z, les quantités X et / ne devront pas varier, 
conformément à la théorie connue des constantes arbitraires ; 

si donc on fait ^=z', et qu’on prenne simplement ~ = m y 
il en résultera 

z’ = m\cos. (y —/). 

On conclut de ces deux équations : 

Xsin. /=-sin. v — z cos. v, 

m 1 

X cos. /=— cos. v -+- z sin. v. 
m 

Donc en appelant 6, et <]>, les termes correspondants de 0' 
et <]/, négligeant les carrés et le produit de z et z\ et mettant 
h et £? à la place de 6 et on aura, d’après les formules du 
n° précédent, 

^ — ~ cos. (y H- C 0 + z sin. (y -f- Çf), 

^ sin. (y +Zt )— z cos.(v -h £*); 

et il ne restera qu’à ajouter ces valeurs avec celles de 0 et Y 
aux formules (20), pour avoir les expressions complètes de 
6 ' et i|/, c’est-à-dire les expressions de ces angles qui répon¬ 
dent au plan tangent à la trajectoire du soleil, mené par le 
centre de la terre. 

Par exemple, la principale inégalité du soleil en latitude, 
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est celle qui provient de l’action de la lune sur la terre. 
En y ayant egard, on a 

z—asm. ( y '—/'), z'—am' cos. (V — /') ; 

a désignant un coefficient constant qui dépend de la masse 
de la lune, et dont la valeur est o ",587 en prenant pour 
cette masse un soixante-quinzième de la terre. On a aussi 

ni =m 13,369; 

au moyen de quoi, l’on trouve 

6 , = 4",a 17 cos. (y 1 — V—v — £ £) -h 3",634 cos. (y 1 — /' 4- v 4- £ £), 
<]/,=—4 ",2 17 sin. ( y '— l — v — ’( t) 4 - 3" ,634 sin.(7/ — /'-h v+Çt). 

( 25 ) D’après le sens suivant lequel les angles <j> et sont 
comptés, le mouvement des équinoxes sera rétrograde, si 
ces angles croissent avec le temps; or, le premier terme de 
la valeur de surpassant les autres, et son coëfficient £ étant 
positif en vertu de l’équation (17), il en résulte que ce mou¬ 
vement est rétrograde, au moins sur lecliptique fixe; et il 
lest effectivement sur les deux écliptiques. Quant à la gran¬ 
deur de ce mouvement, résultant de la valeur de£, elle dé¬ 
pend des moments d’inertie delà terre, et ne peut être dé¬ 
terminée que par l’observation. Pour l’obtenir, nous ferons 
usage de la précession annuelle déterminée par M. Bessel, 
qui l’a trouvée égale à 5 o",37672 pour l’année 170Ô, en la 
rapportant à l’écliptique fixe de la même époque Si l’on 
prend pour unité l’année julienne de 365 J , 25 , et que l’on 
fixe l’origine du temps t au premier janvier 1761 , le coëffi¬ 
cient £ de t dans la valeur de ^ sera cette précession; ainsi 
nous aurons 


ç = 5 o", 37672. 
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Nous emploierons le rapport w de l’action de la lune à celle 
du soleil, que M. Laplace a conclu de l’observation des ma¬ 
rées (*), savoir : 

( 0 = 2 , 35333 . 


Ces valeurs sont les deux données spéciales de la question 
qui nous occupe. On a en outre, selon M. Bessel, 


On a aussi 


h= 23°28'i8". 
e, = o,oi 68 i 4 . 


Quant aux autres constantes contenues dans les formules (18) 
et (20), elles dépendent des masses des planètes que l’on 
adopte. Nous emploierons les valeurs des masses les plus 
récemment déterminées, telles quelles sont données dans la 
cinquième édition de Y Exposition du système du monde. En 
recalculant d’après ces masses, les constantes du tome III 
de la Mécanique céleste , M. Bouvard a trouvé 

f — — o, 000000424696 , 
g , =o“,o 663 i 4 , 

Æ = o", 00001 8658 , 
g'=— o",456917^ 
k ' = o", 00000574 1 . 


Au moyen de ces différentes valeurs, les formules (18) et (20) 
deviennent 


9 = 23 ° 28'18" t\o", 00008001, 
ÿ=t. 5 o", 3 j 5 j 2 — t 2 .o’ jOOO 10906 , 


(*) Mécanique céleste , tome V, page 208. 
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® 2 8'i 8 " —g",45692 — t* .0^000002242 -, 

Y=t. 5 o", 223 oo + v . o",ooo 11637. 

L’angle ô' est l’obliquité moyenne de l’écliptique. Sa gran¬ 
deur pour tannée i 8 i 3 , conclue des solstices d’été de cette 
annee et de 1812 et 1814 observés à Paris, en prenant leur 
résultat moyen et 9 ", 4 o pour le coefficient de la nutation , est 

ô' = 23 ° 27'4 9", 28. 

Si Ion fait t= 6 z ,5 dans l’expression précédente de ô', on 
trouve 

6 = 23 ° 27'49", o 3 . 

La différence entre le calcul et l’observation ne serait donc 
que de o",a 5 ; mais pour vérifier convenablement l’expres¬ 
sion de e', il faudra la comparer à un plus grand nombre de 
solstices consécutifs. Les trois solstices d’hiver observés à la 
même époque donnent une obliquité trop petite de i",45: 
et 1 on sait qu’en effet, ils sont moins propres que ceux d’été 
à cette détermination. 

Si l’on fait t= —2800, on aura l’obliquité qui avait lieu 
1100 ans avant lere chrétienne, ce qui répond à l’époque de 
l’observation chinoise, citée dans l ’Exposition du système du ’ 
monde (*). Suivant cette observation, on aurait 

»’= 23 - 54 ', 

et d après la formule précédente, on a seulement 

ô , =23°49 4 2,r - 


0 Quatrième édition, page 45o. 
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La différence 4 'i 8 " est à peu près le cinquième de la di¬ 
minution totale de l’obliquité depuis cette ancienne époque 
jusqu’à nos jours; elle semblerait indiquer un decroissement 
de l’obliquité plus rapide que nous ne le supposons ; mais il 
est plus probable que cette différence provient de l’inexacti¬ 
tude de l’observation chinoise. 

La précession annuelle est la différence des valeurs de + ' 
qui répondent à t et t -+- i ,* elle est par conséquent 

5 O", 223 ii + £.o"ooo 23274 . 

L’année sidérale est constante et égale en jours moyens à 
365 j , 256374; pendant cet intervalle de temps, le soleil décrit 
36 o°, d’où l’on conclut par une proportion le temps quil 
emploie à décrire l’arc de la précession ; en le retranchant 
ensuite de l’année sidérale, on aura la longueur de l’année 
équinoxiale, savoir : 

365 j ,2422i9 — i.o\ ooooo 6655 ; 

i désignant un nombre de siècles dont 1 origine est toujours 
en 1750. Il en résulte que l’année équinoxiale diminue d a- 
peu près une demi-seconde par siècle. 

(26) Les observations donnent 

— = 0,0748 ^ c=0,08983, € = — 0,33782. 

cette vitesse g du nœud lunaire répondant à l’année julienne 
prise pour unité. La quantité m rapportée à la même unité 
est un peu moindre que 36 o° ; mais 011 peut négliger la dif¬ 
férence dans les formules (19). Au moyen de ces valeurs et de 
celles de Ç, h et w du n° précédent, on trouve 



AUTOUR DE SON CENTRE DE GRAVIT^. 


63 


©=9", 4 oo 4 1. cos. V — o'\oc) 167. cos. 2 V 

-h o", 619 . COS. -h d\oQ 2 . COS. 2 V ’, 

— 17",66677 • sin - *' + o", 84445 . sin. a , 

— 1 ", 196.sin. 21;— o",211 .sin.2 V ; 

V étant la longitude du nœud ascendant de la lune, comptée 
de l’équinoxe du printemps. 

Le coefficient du terme principal de la nutation en latitude 
ou de la valeur de0, est, comme on voit, à très-peu près 9",4o; 
sa valeur, déterminée par l’observation directe, serait 9", 
suivant M. Brinkley, et seulement 8",977 d’après M. de Lin- 
deneau. La différence de ces trois résultats tient en partie aux 
erreurs d.es observations, et aussi à ce que la valeur de <0 que 
nous avons employée, est peut-etre encore un peu trop forte. 
Pour accorder la théorie et le résultat de M. de Lindeneau, 
il faudrait prendre 

w = 2,o3745 ; 

ce qui réduirait la masse de la lune à un 87 e de celle de la 
terre; réduction qu’il est difficile d’admettre, si l’on a égard 
à d’autres phénomènes qui s’accordent à donner une masse 
de la lune plus considérable (*). La valeur de co du n° pré¬ 
cédent suppose cette masse un j 5 e de celle de la terre. 

M. Bessel a le premier introduit dans les calculs astrono¬ 
miques le terme de 0 qui dépend de l’angle 2/'; il n’y a 
pas de raison pour négliger le terme correspondant de W 
dont le coefficient est neuf fois plus considérable. 


(*) Mécanique céleste , tome III, page 159. 


9 
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(27) Il nous reste à considérer la variation du jour solaire 
produite par les déplacements de l’équateur et de l’écliptique. 

Si l’on regarde l’équateur comme fixe, et que l’on trans¬ 
porte au soleil les mouvements de ce plan en sens contraire 
de leurs directions, le jour solaire se réglera sur le mouve¬ 
ment résultant du soleil parallèlement à l’équateur lixe, com¬ 
biné avec la rotation de la terre. Çe mouvement du soleil 
aura lieu dans un pian mobile dont la position sera déter¬ 
minée à chaque instant, par son inclinaison 6' sur le plan 
fixe et par l’angle 8 du n° a 3 , lequel angle exprimera la 
distance du nœud ascendant du soleil sur le même plan fixe, 
à l’intersection de ce plan et de l’écliptique fixe. Désignons par 
a l’ascension droite du soleil, comptée de cette intersection, 
et par v x sa distance angulaire à son .nœud ascendant; nous 
aurons 

tang. (14—£) = cos. 0 ' tang. v x . 

Soit aussi le troisième côté du triangle sphérique du n° cité, 
dont les deux autres côtés sont £ et ^ 4- l, et les angles op¬ 
posés, 6,A et le supplément de 6\ Nous aurons de même 

tang. (v ——/) = cos.'Xtang. (y, — è t ); 

v étant toujours la longitude du soleil comptée sur l’éclip¬ 
tique fixe à partir de la même origine que u. Par les for¬ 
mules connues, on tire de ces deux équations : 

u=$ -h v t — tang.*70'sin.2v I + 7tang. 4 7Ô , sin.4 / v I 

— jtang. 6 jô'sin. 6o> x -t-etc., 
1), = — / H- V + 7 A 2 sill. 2 (y — <[;—/), 

en négligeant les puissances de A supérieures au carré. Le 
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triangle du n° a3 donne, en outre 

Sin .s = “4++i), 

sin. o sin.o 

On a d’ailleurs, d’après ce n°, 

6 =6'—Xsin. -+- /)—|X* sin.*(<|> + /)cot. 6'; 

d’où il résulte 

sin.£,= sin. (<|> h- l) — cos. -f- /)sin. + /) cot. 6' 

— fx’sin. 3 (<|> + /)cot.*6'—7Vcos. (<|/ -4- l) sin.*(<|/ + /); 

et si nous faisons 

Xsin. (4» -+- ï)—p y xcos. (<|> -+- l)=q, 

et que nous négligions toujours le cube de x, nous aurons 
simplement 



$,==<]/ + l—p cot. 0 1 — \pq> 

La valeur de v, deviendra 

v, = v — p cot. 6' — '-p q -+- j (<7* — /?*)sin. 2 v + 7 pq cos. zv ; 
et celle de u prendra la forme : 

u=v — pcot. j9'— '-pq -h P, 
en faisant pour abréger, 

P=7(#*— p') sin.2^ + 7 pqC0&.2.V —tang.^ô'sin.21;, 

-4- 7tang. 4 76'sin. 4^*—f tang. 6 f 6'sin. 6 Vl -h etc. 

La vitesse n de rotation de la terre étant constante, et l’in- 

9 - 
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tégrale jndt , proportionnelle au temps et égalé à nt(n° 1i), 

il en résulte que la distance variable d’un méridien quelcon¬ 
que à une position fixe, sera nt-\-c, en désignant par c 
une constante arbitraire. Si donc on désigne par s l’ascen¬ 
sion droite du soleil à partir de ce méridien, et qu’on fasse 

„ v = m t + e -t- Q, 

m t + s étant la longitude moyenne du soleil, et Q la partie 
périodique de sa longitude vraie, on aura 

s=(n — m) t + p cot. ~ e' + ±pq —P — Q, 

en faisant la constante c —e égale à zéro. Appelons s 'la partie 
non-périodique de cette valeur de s. Pour l’obtenir, il faudra 
faire (n° ^3) : 

p=gt+{k + g'Ut', C 0 t.^’=C 0 t.-,h-^£ Vi , pq=gg' 

et supprimer les termes P et Q, ce qui donnera 
s' = (n — m + g cot. ± 7i) t + Ht 2 , 

OÙ l’on a fait, aussi pour abréger, 

(A +g£—Tgg'cot. j h) cot. l -h = H. 

Le temps apparent et le temps moyen dépendront respec¬ 
tivement des mouvements des astres dont les ascensions droi¬ 
tes seraient ^ et s. La différence de ces deux angles, réduite 
en temps à raison d’une heure pour i5°, exprimera la diffé¬ 
rence du midi vrai au midi moyen, ou ce qu’on appelle 
1 équation du temps. En calculant ses valeurs numériques, 
il faudra employer les valeurs complètes de 6 \ p et q, résul- 
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tant du n° a4 ; mais nous ne nous occuperons ici que de la 
valeur de s' ou du temps moyen. 

(28) Le jour moyen sera l’intervalle de temps pendant 
lequel cet angle s 1 augmentera de 36 o°; en appelant x sa Ion- 
gueur, on aura donc 


36 o°=(rc — m 4 -g-cot. 7 h)x 4- 2x fH. 

L’unité de temps est arbitraire ; si l’on veut qu’elle soit le 
jour moyen à l’époque de i 7 5 o, on aura x= 1 et f=o, et 
par conséquent 

36 o°=rc — m + gcot.{ h. 


La valeur de m donnée par l’observation et rapportée à cette 
unité, est 


m = o°,ÿ 856 i. 


La valeur de g- du n” a 5 répondant à l’année julienne, il fau¬ 
dra la diviser par 365 , 25 , ce qui rendra le terme gcot.,-À 
inutile à considérer. On aura donc 


7*=36 o °, 98661 ; 

et le jour sidéral exprimé en fraction du jour moyen sera 

oi i. 997 2 7 6 * 

Apres avoir divisé les valeurs numériques de g, g et £ 
par 365 ,^ 5 , et celle de k par le carré de ce nombre, afin de 
les rapporter à la nouvelle unité de temps, on trouve 


«_©",0001^344 

(365,2&J* 

Soit en outre 36525 , de manière que i représente un 
nombre de siècles dont 1 origine est en i 7 5o; la grandeur 
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variable du jour moyen sera 

__ *\o,'t328 

X 1 (iooooo) a ’ 

ce qui montre que sa diminution séculaire sera presque in¬ 
sensible. 

Si nous faisons / = t .36 o °, la variable r sera la mesure 
du temps en jours moyens, et d’après les valeurs de s' et de 
H, on aura 

t’.o,oooooi338 
f ~ T (365a5^ 

Le temps t n’est donc pas rigoureusement proportionnel à 
cette mesure; et quoiqu’il s’en écarte fort peu, il sera bon 
cependant d’avoir égard à la différence dans le mouvement 
de la lune,à cause de sa rapidité. Ainsi la longitude moyenne 
m't + z 1 de la lune deviendra 

m 't + e' — i* Tïi (o,oooooi 338 ). 

Le dernier terme s’ajoutera à son équation séculaire, com¬ 
prise dans z ; et le moyen mouvement diurne de la lune étant 

Tnï 13°, i y63, 

à l’époque actuelle, il en résultera 

— i*. o", o 6348 , 

pour la valeur de ce dernier terme ; ce qui diminuera l’équa¬ 
tion séculaire d’environ un 3 oo e de sa valeur. L’effet de cette 
diminution sera, par exemple, d’augmenter d’environ 4o", 
la longitude de la lune en l’an 720 avant notre ère, époque 
des plus anciennes éclipses observées par les Chaldeens, 
qui nous aient été transmises. La même cause produira aussi, 
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dans la longitude du soleil, un très-petit terme proportionnel 
au carré du nombre de siècles, dont le coefficient ne sera 
que d’un millième de seconde. 

Au reste, ces variations séculaires du jour moyen et du 
temps mesuré en jours moyens, comme celles de la préces¬ 
sion annuelle et de l’obliquité de l’écliptique, sont des inéga¬ 
lités périodiques, mais dont la période et l’étendue ne pour¬ 
raient être bien appréciées actuellement, vu l’incertitude qui 
existe encore sur quelques-unes des données de l’observa¬ 
tion dont elles dépendent. Les premiers termes de leurs dé¬ 
veloppements suivant les puissances du temps suffisent aux 
besoins de l’astronomie. 

(29) En substituant dans lequation (17) les valeurs deÇ, 
", h et e < du n ° 25 i observant que la valeur de m rapportée 
h la même unité de temps que celle de £, est 

/72= = 35 9 % 99 3 7 1 ^ 

et que l’on a en outre, 

“=o,oo 2 7 3 o 3 , e'=o,o 548 ü 5 , c=. 0,08983, 
cette équation (17) donne 
2 C —A—B 

-(j-= o,oo65i22. (21) 

Comme on a à tres-peu prèsA = B, il en résulte 
= o,oo 3256 ‘ 1 ; 

et en négligeant le cube de cette fraction, on en déduit 


C—a C —A r . C—A\ 

a c \ 1 H c —) = °^oo3a667. 
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On peut prendre cette valeur pour celle du produit des quan¬ 
tités a et b du n° 2 ; en comparant alors ce produit au rapport 

on voit qu’il le surpasse de plus d’un septième, ainsi 

que nous l’avons dit dans le premier paragraphe de ce Mé¬ 
moire. 

Au moyen des équations (1 4 ) et (21), nous aurons main¬ 
tenant 

= 0,001954. 

L’aplatissement y) de la terre qui paraît résulter de l’ensem¬ 
ble des observations du pendule, est 

7,== 'iï9’ 

ce qui est la valeur du rapport 7 de la force centrifuge à la 
pesanteur sous l’équateur. On aura par conséquent 

et l’on en conclura 

<7 = 1,12922. 

La terre étant un sphéroïde aplati, si elle était en outre 
homogène, le rapport que représente <? (n° i 3 ) serait l’unité ; 
sa masse restant la même, si sa densité est croissante de la 
surface au centre, ce rapport augmentera et pourra devenir 
plus grand que l’unité, comme on le trouve en effet ; ce qui 
montre que le cas d’une densité croissante a réellement lieu. 
D’après l'expérience de Cavendish et les attractions des mon¬ 
tagnes, la densité moyenne de la terre est quintuple de celle 
de l’eau, et à peu près double de celle des substances solides 
qui se rencontrent le plus généralement à sa surface ; or ,on 
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peut satisfaire a cette condition d’une densité superficielle 
moitié de la densité moyenne, et à la valeur précédente de 
g par une infinité d’hypothèses sur la variation de la densité 
intérieure, dont la loi nous est inconnue : je ferai seulement 
remarquer qu’on ne pourrait pas supposer que cette variation 
ne s étendit qua une profondeur peu considérable au-des¬ 
sous de la surface. En effet, si la terre était formée d’un 
noyau sensiblement homogène, recouvert d’une couche dont 
la densité varierait du simple au double dans le sens de son 
épaisseur, en sorte qu’on put la considérer comme une couche 
homogène, d’une densité égale aux trois quarts de la moyenne, 
on trouverait par un calcul facile que pour représenter la va¬ 
leur de <7 dans cette hypothèse, il faudrait que cette épaisseur 
surpassât le quart du rayon total qui est par conséquent 
la moindre profondeur où l’accroissement de densité puisse 
s’étendre. 


io 
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